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Abstract 

We present an overview of recent progress in the phenomenological study of 
neutrino masses, lepton flavor mixing and CP violation. We concentrate on the 
model-independent properties of massive neutrinos, both in vacuum and in matter. 
Current experimental constraints on the neutrino mass spectrum and the lepton 
flavor mixing parameters are summarized. The Dirac- and Majorana-like phases 
of CP violation, which are associated respectively with the long-baseline neutrino 
oscillations and the neutrinoless double beta decay, are discussed in detail. The 
seesaw mechanism, the leptogenesis scenario and the strategies to construct lepton 
mass matrices are briefly described. The features of flavor mixing between one 
sterile neutrino and three active neutrinos are also explored. 
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1 Introduction 



Neutrino physics is a peculiar part of flavor physics, because neutrinos belong to super- 
light flavors in comparison with charged leptons and quarks. There are three central 
concepts in flavor physics: mass, flavor mixing and CP violation. 

• Mass represents the inertial energy possessed by a particle when it exists at rest. A 
massless particle has no way to exist at rest - instead, it must always move at the 
speed of light. A massive fermion (lepton or quark) must exist in both left-handed 
and right-handed states, because the field operators responsible for the nonzero 
mass of a fermion have to be bilinear products of the spinor fields which flip the 
fermion's handedness. 

• Flavor mixing measures the mismatch between flavor eigenstates and mass eigen- 
states of leptons or quarks, caused by the Higgs interactions. Flavor eigenstates 
are the members of the weak isospin doublets that transform into each other via 
the interaction with W ± bosons, while mass eigenstates are the states of definite 
masses created by the interaction with Higgs bosons. If neutrinos were massless, 
lepton flavor mixing would not exist. 

• CP violation means that matter and antimatter or a reaction and its CP-conjugated 
process are distinguishable. It may manifest itself in the weak interactions through 
nontrivial complex phases residing in the flavor mixing matrix of leptons or quarks. 
If CPT is invariant, CP violation will give rise to T violation, or vice versa. 

The problems of fermion masses, flavor mixing and CP violation are correlated with 
one another and fundamentally important in particle physics. The study of these prob- 
lems will ultimately help us understand the nature of matter and the matter-antimatter 
asymmetry of the universe. Since the early 1960's, a lot of experimental efforts have 
been made to measure the parameters of quark mixing and CP violation. The problem 
of lepton mixing and CP violation was not at the front of experimental particle physics 
for a long time, however, partly because neutrinos were assumed to be massless in the 
successful theory of weak and electromagnetic interactions - the standard electroweak 
model [1]. This situation has dramatically changed since the compelling evidence in 
favor of atmospheric and solar neutrino oscillations was achieved by Super-Kamiokande 
(SK) [2], SNO [3], KamLAND [4] and K2K [5] Collaborations in the past five years. 

Neutrino oscillation is a quantum phenomenon which can naturally occur if neutrinos 
are massive and lepton flavors are mixed. Thanks to the elegant SK, SNO, KamLAND 
and K2K experiments, we are now convinced that the long-standing solar neutrino (u e ) 
deficit and the atmospheric neutrino (z/ M ) anomaly are both due to neutrino oscillations. 
The study of neutrino masses and lepton flavor mixing is therefore becoming one of the 
hottest fronts of today's particle physics and cosmology. 

Our present knowledge on the properties of neutrinos comes not only from a number 
of neutrino oscillation experiments, but also from the direct-mass-search experiments, 
the neutrinoless double beta decay experiments and the astrophysical or cosmological 
observations [6]. The following is a partial list of what we have known about neutrino 
masses and lepton flavor mixing. 
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Neutrinos are massive but their masses are tiny. The mass scale of three active 
neutrinos (v e , and u T ) is expected to be at or below 0(1) eV. 



• Two independent mass-squared differences of three neutrinos, which are associated 
separately with solar and atmospheric neutrino oscillations, are very small and have 
a strong hierarchy. Typically, Am^ ~ C(10" 5 ) eV 2 and Am 2 tm ~ 0(1(T 3 ) eV 2 
hold. 

• Two lepton flavor mixing angles, which are related separately to solar and atmo- 
spheric neutrino oscillations, are much larger than the Cabibbo angle of quark 
mixing (0 C « 13°). Typically, 6 sun ~ 33° and atm ~ 45° hold. 

• The other lepton flavor mixing angle, which is relevant to the CHOOZ reactor 
experiment for neutrino oscillations [7], is very small and even vanishing. The 
generous upper bound of this angle is 6* c h z < #c a ^ present. 

It is obvious that the mass spectrum of neutrinos and the mixing pattern of lepton flavors 
are very different from those of quarks. 

However, there exist many open questions about massive neutrinos and lepton flavor 
mixing. For example, 

• Are massive neutrinos Dirac or Majorana particles? If massive neutrinos are Dirac 
particles, they can be distinguished from their antiparticles. By definition, a Majo- 
rana neutrino is identical to its antiparticle. It is possible to identify the Majorana 
nature of massive neutrinos through the observation of the neutrinoless double beta 
decay of some even-even nuclei, in which the total lepton number is not conserved. 

• How many neutrino species are there? We have known that there are three species 
of active neutrinos (u e , and v T ), corresponding to three species of charged lep- 
tons. It remains unclear whether the light sterile neutrinos, which have been 
assumed to interpret the controversial LSND data [8] on I7 M — > V e oscillations, are 
really existent or not. 

• Why are neutrino masses so tiny? The fact that the masses of neutrinos are 
considerably smaller than the masses of charged leptons or quarks is a big puzzle 
to particle physicists. Although a lot of theoretical models about neutrino masses 
have been proposed at either low or high energy scales, none of them has proved 
to be very successful and conceivable. 

• What is the absolute scale of neutrino masses? It is very important to know the 
absolute values of three neutrino masses, because they are fundamental parameters 
of flavor physics. The mass scale of neutrinos is likely to indicate the energy scale 
of new physics responsible for the generation of neutrino masses and lepton flavor 
mixing. Unfortunately, we are only aware of the upper bounds of neutrino masses. 

• Why are the mixing angles 6> sun and 6> atm so big? The bi-large neutrino mixing 
pattern is also a mystery to many theorists, because it is "anomalously" different 
from the familiar tri-small quark mixing pattern. Although the lepton flavor mixing 
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angles are in general expected to relate to the mass spectra of charged leptons and 
neutrinos, their specific relations have not convincingly or model-independently 
been established. 

• How small is the mixing angle # c h z ? Current experimental data yield an upper 
bound 9 chz < 9 C , but the possibility # c h z = 0° cannot be excluded. The smallness 
of 6* c h z requires a good theoretical reason, so does the vanishing of 9 chz . If 9 chz = 0° 
held, there would be no chance to observe leptonic CP or T violation in normal 
neutrino-neutrino and antineutrino-antineutrino oscillations. 

• Is there leptonic CP violation? A necessary condition for the existence of CP and 
T violation in normal neutrino oscillations is 9 chz ^ 0. As CP violation has been 
discovered in the quark sector, we feel that there is no reason why CP should be 
conserved in the lepton sector. It is very difficult to detect the effects of CP or T 
violation in any realistic neutrino oscillation experiments, however. 

• Can the leptonic CP-violating phases be determined? If lepton flavor mixing is 
correlated with CP violation, one has to determine the relevant CP-violating phases 
through various possible experiments. The neutrinoless double beta decay and 
long-baseline appearance neutrino oscillations are expected to be sensitive to the 
Majorana- and Dirac-like phases of CP violation, respectively. To implement such 
measurements remains a big challenge to experimentalists. 

Of course, much more experimental and theoretical efforts are needed to make, in order to 
answer the important questions listed above. Much more phenomenological attempts are 
also needed to make, so as to bridge the gap between experimental data and theoretical 
models. 

The purpose of this article is to give an overview of recent progress in the phenomeno- 
logical study of neutrino masses, lepton flavor mixing and CP violation. We concentrate 
on the model-independent properties of massive neutrinos, although it is unavoidable to 
introduce two very attractive ansatze in today's neutrino physics - the seesaw mechanism 
and the thermal leptogenesis scenario. 

The remaining parts of this article are organized as follows. Section 2 is devoted to a 
review of the neutrino mass spectrum. First of all, we make a short introduction to the 
Dirac and Majorana neutrino masses, and to the well-known seesaw mechanism. Then 
current experimental constraints on neutrino masses, including those from the neutrino- 
less double beta decay, kinematic measurements, neutrino oscillations and cosmological 
observations, are briefly summarized. Finally some comments are given on the strategies 
to construct the phenomenological textures of lepton mass matrices. 

In section 3, the model-independent features of lepton flavor mixing are illustrated 
both in vacuum and in matter. We present a classification of various parametrizations 
of the 3x3 lepton flavor mixing matrix, and highlight one of them which is particularly 
useful for the study of neutrino oscillations. A few concise sum rules for neutrino masses 
and lepton flavor mixing in matter are derived. Several constant mixing patterns of 
massive neutrinos are introduced. Finally we point out the differences and similarities 
between the phenomenon of lepton flavor mixing and that of quark flavor mixing. 
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Section 4 is devoted to leptonic CP and T violation. The rephasing invariants of CP 
violation and the unitarity triangles are described both in vacuum and in matter. A 
few salient features of CP and T violation in long-baseline neutrino oscillations are also 
discussed. We pay some special interest to the thermal leptogenesis scenario to interpret 
the observed baryon-antibaryon asymmetry of the universe. Finally we comment on 
the possible connection or disconnection between the leptonic CP-violating quantities at 
high and low energy scales. 

The conclusion and outlook are presented in section 5. 

It is worth remarking that the main body of this article deals with neutrino masses, 
lepton flavor mixing and CP violation in the scheme of three lepton families. To be 
complete, the properties of flavor mixing between one sterile neutrino (i/ s ) and three 
active neutrinos (v e , and u T ) are discussed in appendix A. 



2 Neutrino Mass Spectrum 
2.1 Theoretical Preliminaries 

The standard model of electromagnetic and weak interactions is based on the gauge 
group SU(2)l x U(1)y [1]- In this framework only the left-handed leptons (and quarks), 
which transform as SU(2) doublets, take part in the charged-current weak interactions: 



e 




-£ cc = ^= (^^Olt" ^ I K + h - c - > (2- 1 ) 

where (e, /z, r) and (z/ e , u T ) are the flavor eigenstates of charged leptons and neutrinos, 
respectively. If SU(2) L were an exact symmetry, all leptons would be massless and their 
flavor eigenstates would be physically indistinguishable from their mass eigenstates. In 
reality, however, this gauge symmetry is badly broken. A spontaneous breakdown of the 
SU(2)l symmetry is realized in the standard model by means of the Higgs mechanism. 
After the symmetry breaking, the charged leptons (and quarks) acquire their masses 
through the Yukawa interactions: 



-d = (e, pi, r) L M x | fi | + h.c. , (2.2) 

R 

where M; denotes the charged lepton mass matrix and its scale is characterized by 
the electroweak symmetry breaking scale v 174 GeV. The vanishing of three neu- 
trino masses follows as a straightforward consequence of the symmetry structure of the 
standard model, in which only a single Higgs doublet exists and the lepton number 
conservation is assumed 1 . 



*It is actually the (B — L) symmetry that makes neutrinos exactly massless in the standard model, 
where B denotes the baryon number and L stands for the total lepton number. The reason is simply 
that a neutrino v and an antineutrino V have different values of (B — L) [9] . Thus the naive argument 
for massless neutrinos is valid to all orders in perturbation and non-perturbation theories, if (B — L) is 
an exact symmetry. 
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Note that the mass eigenstates of charged leptons can always be chosen to coincide 
with their flavor eigenstates through an appropriate but physically-irrelevant unitary 
transformation of the right-handed fields. In such a specific flavor basis the coincidence 
between mass and flavor eigenstates of neutrinos can also be achieved, if neutrinos are 
assumed to be exactly massless Weyl particles. Hence there is no lepton flavor mixing 
within the framework of the standard electroweak model. 



2.1.1 Dirac and Majorana Masses 

However, the assumption of lepton number conservation or masslessness of neutrinos 
is not assured by any basic symmetry principle of particle physics. Most reasonable 
extensions of the standard model (such as the grand unified theories) do allow the absence 
of lepton number conservation and the existence of nonvanishing neutrino masses. If 
neutrinos are really massive and their masses are non-degenerate, it will in general be 
impossible to find a basis of the flavor space in which the coincidence between flavor 
and mass eigenstates holds both for charged leptons and for neutrinos. In other words, 
the flavor mixing phenomenon is naturally expected to appear between three charged 
leptons and three massive neutrinos, just like the flavor mixing between three up-type 
quarks and three down- type quarks [10]. 

If neutrinos have nonzero masses, they may be either Dirac or Majorana particles. 
The field of a massive Dirac neutrino describes four independent states: left-handed 
and right-handed particle states (u L and v R ) as well as left-handed and right-handed 
antiparticle states (P L and P R ). Among them the u L and P R states, which already exist 
in the standard model, can take part in weak interactions. The v R and P L states need 
to be introduced into the standard model as necessary ingredients to give the Dirac 
neutrino a mass, but they should be "sterile" in the sense that they do not take part 
in the normal weak interactions. A Dirac mass term, which conserves the total lepton 
number but violates the law of individual lepton flavor conservation, can be written as 
follows: 

/ 

v» I + h.c. , (2.3) 



-£ D = {ye, Vp, V T ) L M D 

V^/r 

where M D denotes the 3x3 Dirac neutrino mass matrix. The mass term in (2.3) is 
quite similar to the mass term of charged leptons in (2.2), hence the scale of Md should 
also be characterized by the gauge symmetry breaking scale v. In this case the Yukawa 
coupling constants of three neutrinos must be extremely smaller than those of three 
charged leptons or six quarks, such that tiny neutrino masses can result. This dramatic 
difference between the Yukawa couplings of neutrinos and charged leptons (or quarks) is 
commonly considered to be very unnatural in a sound theory of fermion mass generation, 
however. 

On the other hand, the neutrino v may be a Majorana particle, which has only 
two independent states of the same mass (u L and P R , or u R and P L ). By definition, a 
Majorana neutrino is its own antiparticle: v c = Cu T = e l ®v [11], where C denotes the 
charge-conjugation operator and 6 is an arbitrary real phase. A Majorana mass term, 
which violates both the law of total lepton number conservation and that of individual 



7 



lepton flavor conservation, can be written either as 



-£m ( l) = g v v» v t)l M l I vl ] + h.c. , (2.4) 



or as 



T / R 



-^M(R) = ^ (^,^,^) L Mr [ I + h.c. , (2.5) 



R 

where M L and M R denote the symmetric 3x3 mass matrices of left-handed and right- 
handed Majorana neutrinos, respectively. Note that the mass term £m(l) cannot natu- 
rally arise from a simple theory of electroweak interactions which is invariant under the 
SU(2)l x U(l)y gauge transformation and has no SU(2)l triplet field (like the case in the 
standard model). It is possible to incorporate the mass term £m(r) into an electroweak 
theory with SU(2) L x U(l)y gauge symmetry, because right-handed neutrinos are SU(2) L 
singlets. In general, a neutrino mass Lagrangian may include all of the above-mentioned 
terms: C^, £m(l) and £m(r), i n which M D , M L and M R are complex mass matrices. 

A variety of specific theoretical and phenomenological models have been prescribed in 
the literature to interpret how lepton masses are generated and why neutrino masses are 
so tiny. Regardless of the energy scales at which those models are built, the mechanisms 
responsible for fermion mass generation and flavor mixing can roughly be classified into 
five different categories [12]: (a) Radiative mechanisms [13]; (b) Texture zeros [14]; (c) 
Family symmetries [15]; (d) Seesaw mechanisms [16]; and (e) Extra dimensions [17]. 
Among them, the seesaw mechanisms are particularly natural and interesting. A brief 
introduction to the seesaw idea will be presented in section 2.1.2. For recent reviews of 
other interesting models and ansatze on neutrino masses, we refer the reader to Ref. [12] 
andRefs. [18]-[21]. 

2.1.2 The Seesaw Mechanism 

A simple extension of the standard model is to include one right-handed neutrino in each 
of three lepton families, while the Lagrangian of electroweak interactions keeps invariant 
under the SU(2) L x U(l)y gauge transformation. After spontaneous symmetry breaking, 
the resultant lepton mass term £ mass consists of £/ (charged lepton term), £ D (Dirac 
neutrino term) and £m(r) (righ-handed Majorana neutrino term), which have been given 
in (2.2), (2.3) and (2.5). To be explicit, we have 

1- 



" £ - = - C ' + 2 ( "'^l«S JMkJUV (2 ' 6) 

where v denotes the column vector of v e , and u T fields; i.e., v T = (v e , v T ) or 
u c = (z/g, 1/° v c T ) T . In obtaining (2.6), we have used the relation V^M^v^ = (c c )lM^(/) r 
as well as the properties of v and z/°. As pointed out in section 2.1.1, the scale of Md 
is characterized by the electroweak symmetry breaking scale v. The scale of Mr can 
naturally be much higher than v, because right-handed neutrinos are SU(2)l singlets and 
their corresponding mass term is not subject to the scale of gauge symmetry breaking. It 
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Figure 2.1: Illustrative plot for the seesaw mass term. 



has commonly been expected that the scale of Mr is not far away from the grand unified 
theory (GUT) scale Aqtjt ~ 10 16 GeV. In this case, the smallness of left-handed neutrino 
masses (m 4 <C v for i — 1, 2, 3) is essentially attributed to the largeness of right-handed 
neutrino masses (Mj 3> t> for i = 1,2,3). Such an elegant idea, the so-called seesaw 
mechanism, was first proposed by Yanagida and by Gell-Mann, Ramond and Slansky in 
1979 [16]. Fig. 2.1 illustrates the seesaw mass term given in (2.6). 

To show how the seesaw mechanism works, we diagonalize the symmetric 6x6 
neutrino mass matrix in (2.6) by use of the following unitary transformation: 



V RY f M D \ /V RY _(M U 
S Uj Im d t M R J [S Uj - { M R 



(2.7) 



where R, S, U and V are the 3x3 matrices; M v and M R denote the 3x3 diagonal mass 
matrices with eigenvalues and Mj (for i = 1,2,3), respectively. It is straightforward 
to obtain 

VM V V T = A y - M D M£ J <, 

UM K U T = Au + M R , (2.8) 

where Ay = V & M£R*U T M^UR) M b S*V t and = M^R*U T + URlM D . Note 
that both R and S in (2.7) are expected to be of order Md/Mr, as a result of the 
huge hierarchy between the scales of Md and Mr. Hence Ay and A^ fa arc 
excellent approximations for (2.8). We then arrive at Mr UM r U t for the heavy 
(right-handed) neutrino mass matrix and the well-known seesaw formula for the light 
(left-handed) neutrino mass matrix 2 , 

M v = VM U V T ps -MdMr 1 ^ . (2.9) 

The seesaw mechanism was actually motivated by the SO (10) GUTs [22], which include 
right-handed neutrinos automatically together with other charged fermions in irreducible 
16-dimensional multiplets. In such an attractive framework, Md can simply be related to 
the up- type quark mass matrix M u , whose possible textures have well been restricted by 
current experimental data [12]. If both Md (= M u ) and Mr were assumed to be diagonal, 
M v would also be diagonal and its three mass eigenvalues would be m\ = m^/Mi, 
m 2 — 777^/M 2 and m 3 = mf/M 3 . In this illustrative case, m 3 ~ 0.1 eV may simply result 
from M 3 ~ 10 14 GeV, a scale not far away from the GUT scale Aqut- 

2 Note that the minus sign on the right-hand side of (2.9) can always be absorbed via a redefinition 
of the approximately unitary matrix V. For instance, VM V V T w MdM^M^ with V = iV holds. 
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Without loss of generality, both the charged lepton mass matrix M/ and the heavy 
Majorana neutrino mass matrix Mr can be arranged to be diagonal and have positive 
mass eigenvalues. In this specific flavor basis, the approximately unitary matrix V in 
(2.9) links the mass eigenstates of light neutrinos to their flavor eigenstates; i.e., it is 
relevant to the phenomenon of lepton flavor mixing at low energies (see section 3.1 for 
more detailed discussions). A full parametrization of V requires three mixing angles and 
three CP-violating phases. Thus M v consists of nine free parameters. Note that the 



Dirac neutrino mass matrix Md can be parametrized as Md = iV\J M v O\/M-& [23], 
in which O is a complex orthogonal matrix containing three rotation angles and three 
phase angles. It seems that Md in this parametrization consists of eighteen parameters: 
six from V, three from M u , six from O and three from Mr. However, only fifteen of 
them (nine real parameters and six CP- violating phases) are independent, due to the 
seesaw relation. We conclude that the lepton mass matrices Mi, Md and Mr totally 
have twenty-one free parameters: three from M;, fifteen from Md and three from Mr in 
the chosen flavor basis 3 . 

It is worth mentioning that (2.9) has been referred to as the Type-I seesaw formula. 
A somehow similar relation, the so-called Type-II seesaw formula [25], can be derived 
from the generalized lepton mass term 

1- 



where Ml has been defined in (2.4). As Ml must result from a new Yukawa-interaction 
term which violates the SU(2) L x U(l)y gauge symmetry in the standard model, its 
scale might be much lower than the electroweak symmetry breaking scale v. The strong 
hierarchy between the scales of Mr and Ml or Md allow us to make some safe approx- 
imations in diagonalizing the 6x6 neutrino mass matrix in (2.10). We find that the 
heavy (right-handed) neutrino mass matrix remains to take the form Mr ps UM^U t . 
In contrast, the light (left-handed) neutrino mass matrix is given by 

M v = VM U V T ss M L - M D M^M^ . (2.11) 

This result is just the Type-II seesaw formula. As emphasized in Ref. [25], the Type-II 
seesaw relation is a reflection of the left-right symmetry of the theory at high energies. For 
the phenomenological study of neutrino masses, however, the Type-I seesaw mechanism 
is more popular and useful because it involves fewer free parameters than the Type-II 
seesaw mechanism. 

Of course, there are some other versions of the seesaw mechanism, which are either 
more complicated [25] or simpler [26] than what we have discussed above. While the 
seesaw idea is qualitatively elegant to interpret the smallness of left-handed neutrino 
masses, it cannot lead to quantitative predictions unless specific textures of neutrino 
mass matrices (e.g., Md and Mr) are assumed. Many interesting ansatze of lepton mass 
matrices have so far been proposed and incorporated with the seesaw mechanism (for 
recent reviews with extensive references, see Ref. [12] and Refs. [18]— [21]), but they are 
strongly model-dependent and thus beyond the scope of the present review article. 

3 The parameter counting is identical in models with and without supersymmetry [24]. 
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2.2 Experimental Constraints 



2.2.1 Neutrino Oscillations 

Neutrino oscillation is a quantum phenomenon which can naturally happen if neutrinos 
are massive and lepton flavors are mixed. In a simple two-neutrino mixing scheme, the 
neutrino flavor eigenstates v a and are linear combinations of the neutrino mass eigen- 
states u a and u b [27] : v a = u a cos 6 + v h sin 9 and z/g = u b cos 9 — u a sin 9, where 9 denotes 
the flavor mixing angle. Then the probabilities of neutrino oscillations are governed 
by two characteristic parameters: one of them is the neutrino mass-squared difference 
Am 2 = m 2 — m 2 and the other is the flavor mixing factor sin 2 29. Corresponding to the 
"disappearance" and "appearance" neutrino experiments, the survival and conversion 
probabilities of a neutrino flavor eigenstate v a can explicitly be expressed as 

P(u a - v a ) = 1 sin 2 20sin 2 (l.27- A? 




P(y a ^v) = sin 2 2# sin 2 1.27— — , (2.12) 



where f3 ^ a, E is the neutrino beam energy (in unit of GeV), and L denotes the 
distance between the neutrino source and the neutrino detector (in unit of km). A 
lot of experimental data, including those from solar, atmospheric and reactor neutrino 
oscillation experiments, have been analyzed by use of (2.12). 

(1) The first model-independent evidence for neutrino oscillations was obtained from 
the Super-Kamiokande (SK) experiment [2] on atmospheric neutrinos, which are pro- 
duced in the earth's atmosphere by cosmic rays and are detected in an underground 
detector. If there were no neutrino oscillation, the atmospheric neutrinos entering and ex- 
iting the detector should have a spherical symmetry. In other words, the downward-going 
and upward-going neutrino fluxes should be equal to each other: § a (9 z ) = & a (n — 9 Z ) 
versus the zenith angle 9 Z (for a = e or /i). The SK Collaboration has observed an 
approximate up-down flux symmetry for the atmospheric v e neutrinos and a significant 
up-down flux asymmetry for the atmospheric neutrinos. For instance, 

«,(-! < «»■ < -0.2) ^ . 54±0 , 04 ^ ! (2 . 13) 
$ M (+0.2 < cos0 z < +1) T v ; 

was measured for the multi-GeV neutrinos, and it is in conflict with the up-down flux 
symmetry. This result can well be interpreted in the assumption of — > u T neutrino 
oscillations. Current SK [2] and CHOOZ [7] data have ruled out the possibility that 
the atmospheric neutrino anomaly is dominantly attributed to — > v e or — > v s 
oscillations, where v s stands for a "sterile" neutrino which does not take part in the 
normal electroweak interactions (see appendix A for detailed discussions). 

The hypothesis of atmospheric neutrino oscillations has recently received very strong 
support from the K2K [5] long-baseline neutrino experiment, in which the beam is 
produced at the KEK accelerator and measured 250 km away at the SK detector. The 
K2K Collaboration observed a reduction of the flux and a distortion of the energy 



11 



spectrum, which must take place in the presence of — > v T oscillations. The possibility 
that the K2K result is due to a statistical fluctuation instead of neutrino oscillations is 
found to be less than 1% [5]. 

It has been shown that the standard z/ M — > u T oscillations provide the best description 
of the combined SK and K2K data [28], from which Am^ tm = (2.6 ± 0.4) x 10~ 3 eV 2 
and sin 2 26 , atm = l.OOl^os are determined at the la level. At the 90% confidence level, 
we have 

1.65 x 10~ 3 eV 2 < Am 2 tm < 3.25 x 10~ 3 eV 2 (2.14) 

and 0.88 < sin 2 20 atm < 1.00. 

(2) The long-standing problem associated with solar neutrinos is that the flux of solar 
v e neutrinos measured in all experiments (e.g., the SK [2] and SNO [3] experiments as 
well as the earlier Homestake [29], GALLEX-GNO [30] and SAGE [31] experiments) is 
significantly smaller than that predicted by the standard solar models [32]. The deficit 
of solar z/ e 's is not the same in different experiments, implying that the relevant physical 
effect is energy-dependent. It has been hypothesized that the solar neutrino problem 
is due to the conversion of solar z/ e 's into other active or sterile neutrinos during their 
travel from the core of the sun to the detectors on the earth. The SNO experiment 
has model-independently demonstrated that v e — > and (or) u e — > v T transitions are 
dominantly responsible for the solar neutrino deficit. 

What the SNO Collaboration has measured is the flux of solar 8 B neutrinos via 
the charged-current (CC), neutral-current (NC) and elastic-scattering (ES) reactions: 
i> e + d ^ e~ + p + p, Ua + d—^Ua+p + n and u a + e~ — > u a + e~, where a = e, fi or r. 
In the presence of flavor conversion, the observed neutrino fluxes in different reactions 
satisfy 

$ cc = $£ $ nc = $£ + $ ) $ es = $g + ^ w $£ + Q.154V , (2.15) 

where 0"^/<r e ~ 0.154 is the ratio of elastic is e -e and v^-pu scattering cross-sections, and 
denotes the flux of active non-electron neutrinos. Of course, = or equivalently 
«£CC _ ^nc _ ^es wou \d holci, if there were no flavor conversion. The SNO data [3] 
yield 

$ e = (1.76 ± 0.10) x 10 6 cm" 2 s" 1 , $ MT = (3.41+°;^) x 10 6 cm- 2 s" 1 , (2.16) 

a convincing evidence (at the 5.3<r level) for the existence of and v T neutrinos in the 
flux of solar neutrinos onto the earth. 

The flavor conversion of solar v e neutrinos is most likely due to neutrino oscillations. 
In the scheme of two-neutrino oscillations, a global analysis of all available experimental 
data on solar neutrinos (in particular, those from the SK and SNO measurements) leads 
to several regions allowed for the parameters Am 2 un and tan 2 sun : the SMA (small mixing 
angle), LMA (large mixing angle) and LOW (low mass-squared difference) regions based 
on the Mikheev-Smirnov-Wolfenstein (MSW) mechanism [33] as well as the VO (vacuum 
oscillation) and other possible regions [34], among which the LMA region is most favored. 
Recently the LMA solution has been singled out as the only acceptable solution to the 
solar neutrino problem, thanks to the KamLAND [4] experiment. 
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The KamLAND Collaboration has measured the flux of F e 's from distant nuclear 
reactor via the inverse /3-decay reaction V e + p — > e + + n. The typical baseline of this 
experiment is 180 km, allowing for a terrestrial test of the LMA solution to the solar 
neutrino problem. The ratio of the observed inverse /5-decay events to the expected 
number without V e disappearance is 0.611 ± 0.094 for E u > 3.4 MeV. Such a deficit can 
naturally be interpreted in the hypothesis of neutrino oscillations, and the parameter 
space of (Am 2 , sin 2 29) is found to be compatible very well with the LMA region. We 
are then led to the conclusion that the LMA solution is the only correct solution to the 
solar neutrino problem. A global analysis of the combined SK, SNO and KamLAND 
data yields [35] 

5.9 x 10~ 5 eV 2 < Am 2 un < 8.8 x 10~ 5 eV 2 (2.17) 

and 0.25 < sin 2 9 sun < 0.40 at the 90% confidence level. 

(3) The purpose of the CHOOZ and Palo Verde reactor experiments [7] is to search 
for V e — > V e disappearance neutrino oscillations in the atmospheric range of Am 2 . No 
indication in favor of neutrino oscillations was found from both experiments, leading to 
a strong constraint on the flavor mixing factor: sin 2 2# chz < 0.10 for Am 2 hz > 3.5 x 
10 -3 eV 2 ; or sin 2 29 chz < 0.18 for Am 2 hz > 2.0 x 10~ 3 eV 2 . The impact of this constraint 
on the lepton flavor mixing matrix will be discussed in section 3.1. 

Let us denote the mass eigenstates of v e , and v T neutrinos as (^1,^2,^3), whose 
eigenvalues are (7711,7712,7713). There are only two independent neutrino mass-squared 
differences: Am^ = m 2 — m 2 and Am 2 2 = m 2 — m 2 or Am^ = m§ — m 2 . Without loss 
of generality, we make the identification 

Am s 2 un = \Am 2 21 \ « |Am 2 2 | = Am 2 tm , (2.18) 

where the big hierarchy between Am 2 un and Am 2 tm has been taken into account. In- 
deed, Am 2 un <C Am 2 tm together with sin 2 29 c \ a <C 1 implies that solar and atmospheric 
neutrino oscillations are approximately decoupled. The deficits of solar and atmospheric 
neutrinos are dominated respectively by v e — > and — > v T transitions. With the 
help of (2.18), mi and m 2 can be given in terms of m 3 , Am 2 un and Am 2 tm [36]: 

mi = ^ml+pAml tm + q Aml Q , 

m 2 = ^mj+pAml tm , (2.19) 



= ±1 stand for four possible neutrino mass spectra : 



where p — ±1 and q 



(p,g) = (-l,-l) 
(p,g) = (-l,+l) 
(p,g) = (+l,-l) 
(p,q) = (+1,+1) 



mi < m 2 < m 3 , 

m l > m 2 < m3 , 
m l < m 2 > m 3 , 

mi > m 2 > 777.3 • 



(2.20) 



Current solar neutrino oscillation data favor q — — 1; i.e., mi < m 2 . The sign of p may 
be determined from the future long-baseline neutrino oscillation experiments. 

4 The mass convention m' x < m! 2 < m' 3 has been used in some literature for the redefined neutrino 
mass eigenstates v[, which are related to the standard neutrino mass eigenstates Vi by an orthogonal 
transformation [37]. We shall not adopt v[ in the present article. 
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Figure 2.2: Illustrative correlation of three neutrino masses for p — ±1 cases. 



Given the best-fit values Am s 2 UI1 = 7.3 x 1(T 5 eV 2 and Am 2 tm = 2.5 x 1(T 3 eV 2 , the 
numerical correlation of m 1 , m 2 and m 3 is illustrated in Fig. 2.2 by use of (2.19). We see 
that three masses become nearly degenerate, if one of them is larger than 0.1 eV (i.e., 
rrii > 0.1 eV). For the p — — 1 case, the lower bound of m 3 is m 3 > ^ Am^ + Am 2 UI1 , 

which leads to m 2 > Am 2 un . A normal hierarchy m x < m 2 < m 3 may appear, if 
mi ~ (9(10~ 2 ) eV or smaller. For the p — +1 case, however, mi « m 2 > m 3 always 
holds, no matter how small m 3 is taken (inverted hierarchy). 

Finally it is worthwhile to mention the LSND evidence [8] for — > v e neutrino 
oscillations with Am1 SND ~ 1 eV 2 and sin 2 2#lsnd ~ 10~ 3 — 10~ 2 . This observation 
was not confirmed by the KARMEN experiment [38], which is sensitive to most of the 
LSND parameter space. The disagreement between these two measurements will be 
resolved by the MiniBOONE experiment [39] at Fermilab in the coming years. Before 
the LSND result is confirmed on a solid ground, the conservative approach is to set it 
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aside tentatively and to concentrate on solar and atmospheric neutrino oscillations in the 
mixing scheme of three lepton families. A natural way to simultaneously accommodate 
solar, atmospheric and LSND data is to assume the existence of a light sterile neutrino, 
which may oscillate with three active neutrinos. Such a four-neutrino mixing scheme 
will be discussed in appendix A. 

2.2.2 Cosmological Bounds 

Neutrinos play an important role in cosmology and astrophysics [40]. Useful information 
on neutrino masses can be obtained, for example, from studies of the cosmological relic 
density and dark matter, the power spectrum in large scale structure surveys, the Lyman 
a forest, the Big Bang nucleosynthesis, the supernovae, the ultra high energy (UHE) 
cosmic rays, and the gamma ray bursts. For recent reviews of these interesting topics, we 
refer the reader to Refs. [40]-[43]. Here we focus our attention on the latest cosmological 
bound on the sum of light neutrino masses, obtained from the impressive data of the 
Wilkinson Microwave Anisotropy Probe (WMAP) [44]. 

According to the Big Bang cosmology, neutrinos were in thermal equilibrium with 
photons, electrons and positrons in the early universe. When the universe cooled to 
temperatures of 0(1) MeV, neutrinos decoupled from the primordial e ± 7 plasma, leading 
to the fact that the present-day number density of neutrinos is similar to that of photons. 
If neutrinos are massive, they contribute to the cosmological matter density [45], 

n uh 2 = V^^, (2.2i) 

Y 93-5 eV ' K J 

where Q u denotes the neutrino mass density relative to the critical energy density of 
the universe, and h is the Hubble constant in units of 100 km/s/Mpc. A tight limit 
Vt u h 2 < 7.6 x 10~ 3 (at 95% C.L.) has recently been extracted from the striking WMAP 
data [44], combined with the additional cosmic microwave background (CMB) data sets 
(CBI and ACBAR) and the observation of large scale structure from 2dF Galaxy Redshift 
Survey (2dFGRS) [46]. This impressive bound leads straightforwardly to 5 

^m, < 0.71 eV , (2.22) 

i 

which holds at the same confidence level. Some discussions are in order. 

(a) A tight upper bound on can be achieved from (2.22) together with current 
data on solar and atmospheric neutrino oscillations. To see this point more clearly, we 
calculate the dependence of mi +m 2 + m 3 on m 3 with the help of (2.19) and (2.22). The 
numerical results are shown in Fig. 2.3, where the best-fit values Am s 2 UI1 = 7.3 x 10~ 5 eV 2 
and Am 2 tm = 2.5 x 10 -3 eV 2 have typically been input. Note that only the q = — 1 or 
mi < m 2 case, which is supported by current solar neutrino oscillation data, is taken 
into account. For the p = — 1 or m 2 < m 3 case, m 3 has an interesting lower bound 

5 The exact value of this upper bound depends on other cosmological parameters, as emphasized by 
Hannestad in Ref. [47], where a somewhat more generous limit mi + m 2 + m 3 < 1.0 eV (at 95% C.L.) 
has been obtained for N v = 3. 



15 



10 



mi + rri2 + m 3 (eV) 
WMAP upper bound 




m 2 > m 3 ^^-^ / 
m 2 < m 3 

i i 





10~ 3 10~ 2 10 _1 i 



m 3 (eV) 

Figure 2.3: Illustrative dependence of mi + m 2 + m 3 on m 3 . The WMAP result sets an 
upper limit on m 3 ; i.e., m 3 < 0.24 eV for both m 3 > m 2 and m 3 < m 2 cases. 



^3 > Y^ m atm + Am s 2 UI1 m 0.051 eV; but for the p = +1 or m2 > m 3 case, even m 3 = 
is allowed (inverted hierarchy). We see that these two cases become indistinguishable for 
m 3 > 0.2 eV, implying the near degeneracy of three neutrino masses. Once the WMAP 
limit in (2.22) is included, we immediately get m 3 < 0.24 eV. As a consequence, we have 
rrii < 0.24 eV for % — 1, 2 or 3. Similar results have been obtained in Refs. [47, 48]. 

(b) The LSND [8] evidence for neutrino oscillations is strongly disfavored, because 
^ m LSND ~ 1 e V 2 is essentially incompatible with (2.22). However, a marginal agreement 
between the LSND data and the WMAP data is not impossible in the four-neutrino 
mixing scheme, if the latter's limit on neutrino masses is loosened to some extent (e.g., 
mi + m 2 + m 3 + m 4 < 1.4 eV for N v = 4 [47]). In this sense, it might be premature to 
claim that the LSND result has been ruled out by the WMAP data [48]. We hope that 
the MiniBOONE experiment [39] can ultimately confirm or disprove the LSND result. 

2.2.3 Kinematic Measurements 

An essentially model-independent way to determine or constrain neutrino masses is to 
measure some typical weak decays whose final-state particles include neutrinos. The 
kinematics of such a process in the case of non-zero neutrino masses is different from 
that in the case of zero neutrino masses. This provides an experimental opportunity to 
probe neutrino masses in a direct way. In practice, direct neutrino mass measurements 
are based on the analysis of the kinematics of charged particles produced together with 
the neutrino flavor eigenstates \u a ) (for a = e,fi,r), which are superpositions of the 
neutrino mass eigenstates |z/j) (for i = 1,2,3): \u a ) = V a i\vi) + V a 2\v 2 ) + V a3 \u 3 ) with V 
being the 3x3 lepton flavor mixing matrix. The effective masses of three neutrinos can 
then be defined as 




(2.23) 



16 



for a = e,fi and r. The kinematic limits on (m) e , (m)^ and (m) T can respectively be 
obtained from the tritium beta decay fH — ^He + e~ +v e , the n + — > /i + + i/^ decay and 
the r — > 57r + z/ r (or r — > 37T + u T ) decay [6]: 

(m) e < 2.2 eV , 
(m) M < 0.19 MeV , 

(m) T < 18.2 MeV . (2.24) 

We see that the experimental sensitivity for (m)^ is more than four orders of magnitude 
smaller than that for (m) e , and the experimental sensitivity for (m) T is two orders of 
magnitude lower than that for (m)^. 

Note that the unitarity of V leads straightforwardly to a simple sum rule between 
{m) 2 a and mf [49]: 

E<™>« = £™? < feWf • (2-25) 

As the sum of three relativistic neutrino masses has well been constrained by the recent 
WMAP data [44], m 1 + m 2 + m ? ,< 0.71 eV at the 95% confidence level, we are led to 

(m)l + (m)l + (m) 2 T < 0.50 eV 2 . (2.26) 

This generous upper bound implies that (m)"^ < 0.50 eV 2 or (m) a < 0.71 eV holds 
for a = e,fi and r. One can see that the cosmological upper bound of (m) e is about 
three times smaller than its kinematic upper bound given in (2.24). The former may 
be accessible in the future KATRIN experiment [50], whose sensitivity is expected to 
be about 0.3 eV. In contrast, the cosmological upper bound of (m)^ is more than five 
orders of magnitude smaller than its kinematic upper bound given in (2.24), and the 
upper limit of (m) T set by the WMAP data is more than seven orders of magnitude 
smaller than its kinematic upper limit. It seems hopeless to improve the sensitivity of 
the kinematic (m)^ and (m) T measurements to the level of 0.71 eV. 

If current neutrino oscillation data are taken into account, however, more stringent 
upper limits can be obtained for the effective neutrino masses (m) e , and (m) T . 

Substituting (2.19) into (2.23), we get 

(m) a = yfml +p(l- |K 3 | 2 ) AmL + q\V al \ 2 Am^ . (2.27) 

The present experimental data indicate \V e z\ 2 ~ sin 2 # cnz 1 [7] and Am 2 un Am 2 tm . 
Furthermore, sin 2 2^ atm ps 1 is strongly favored, leading to jV^I ~ sin^ atm ~ l/y/2 and 
\V r3 \ w cos^ a t m w 1 / y/2 for \V e3 \ < 1. Therefore, (2.27) can be simplified as 

(m) e w \jm\ + pAml tm , 

(™)lM ~ \j m l + | Am atm , 

(m) T » ^jml + |AmL n . (2.28) 

We can see that (m) M « (m) T is a consequence of current neutrino oscillation data. In 
addition, (2.28) shows that (m) e is slightly larger than (m)^ and (m) T for p = +1 or 
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Figure 2.4: Illustrative dependence of (m) a (for a = e,fi,r) on m 3 . Curve a: (m) e with 
m 3 < m 2 ; Curve b: (m)^ ~ (m) T with m 3 < m 2 ; Curve c: (m)^ ~ (m) T with f«3 > m 2 ; 
and Curve d: (m) e with m 3 > m 2 . The WMAP result leads to (m) a < 0.24 eV. 

m 2 > m 3 ; and it is slightly smaller than (m) M and (m) T for p = — 1 or m 2 < m 3 . The 
numerical dependence of (m) e , and (m) T on m 3 is illustrated in Fig. 2.4. We find 
that it is impossible to distinguish between the m 2 < m 3 case and the m 2 > m 3 case for 
m 3 > 0.2 eV. Note that the dependence of {m) a on m 3 is similar to the dependence of 
m \ + m 2 + 1^3 on m 3 shown in Fig. 2.3. Taking account of the upper limit m 3 < 0.24 
eV obtained in section 2.2.2, we arrive at (m) a < 0.24 eV for a = e,fi or r. This upper 
bound is suppressed by a factor of three, compared to the upper bound obtained from 
(2.26) which is independent of the neutrino oscillation data. 

The result (m)^ ps (m) T < 0.24 eV implies that there is no hope to kinematically 
detect the effective masses of muon and tau neutrinos [49] . As the WMAP upper bound 
is in general valid for a sum of the masses of all possible relativistic neutrinos (no matter 
whether they are active or sterile), it seems unlikely to loosen the upper limit of {m) a 
obtained above in the assumption of only active neutrinos. Therefore, the kinematic 
measurements of (m) M and (m) T have little chance to reveal the existence of any exotic 
neutral particles with masses much larger than the light neutrino masses. 

2.2.4 Neutrinoless Double-/? Decay 

The neutrinoless double beta decay of some even-even nuclei, 

A(Z,N) -> A(Z + 2,N-2) + 2e~ , (2.29) 

can occur through the exchange of a Majorana neutrino between two decaying neutrons 
inside a nucleus, as illustrated in Fig. 2.5. It would be forbidden, however, if neutrinos 
were Dirac particles. Thus the neutrinoless double-/? decay provides us with a unique 
opportunity to identify the Majorana nature of massive neutrinos. 
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A(Z,N)^^= Nuclear Physics A(Z+2,N-2) 



Figure 2.5: Illustrative plot for the neutrinoless double beta decay of some even-even 
nuclei via the exchange of a virtual Majorana neutrino. 

The rate of the neutrinoless double-/? decay depends both on an effective neutrino 
mass term (m) ee and on the associated nuclear matrix element. The latter can be 
calculated, but some uncertainties are involved in the calculations [43]. In the three- 
neutrino mixing scheme, the effective Majorana mass (m) ee is given by 

(m) ee = \m 1 V e \+m 2 V^ + m 3 V^\ , (2.30) 

where m ; (for i = 1,2,3) denote the physical masses of three neutrinos, and V ei stand 
for the elements in the first row of the 3x3 lepton flavor mixing matrix V. It is obvious 
that (m)ee = would hold, if mi — were taken. 

Many experiments have been done to search for the neutrinoless double-/? decay [51]. 
Among them, the Heidelberg- Moscow [52] and IGEX [53] 76 Ge experiments have the 
highest sensitivity and yield the most stringent upper bounds on (m) ee : 

(m) e e < (0.35 - 1.24) eV (Heidelberg - Moscow) , 

(m) ee < (0.33 - 1.35) eV (IGEX) , (2.31) 

where the uncertainties associated with the nuclear matrix elements have been taken 
into account 6 . A number of new experiments [51], which may probe (m) ee at the level 
of 10 meV to 50 meV, are in preparation. 

While (ra) ee ^ must imply that neutrinos are Majorana particles, (m) ee = does 
not necessarily imply that neutrinos are Dirac particles. The reason is simply that the 
Majorana phases hidden in V e j may lead to significant cancellations on the right-hand 
side of (2.30), making (m) ee vanishing or too small to be detectable [57]. Hence much 
care has to be taken, if no convincing signal of the neutrinoless double beta decay can 
be experimentally established: it may imply that (1) the experimental sensitivity is 
not high enough; (2) the massive neutrinos are Dirac particles; or (3) the vanishing or 
suppression of (m) ee is due to large cancellations induced by the Majorana CP-violating 
phases. The third possibility has recently been examined [58] from a model-independent 
point of view and with the help of the latest experimental data. It is found that current 
neutrino oscillation data do allow (m) ee = to hold, if the Majorana phases lie in two 

6 After re-analyzing the data from the Heidelberg-Moscow experiment, Klapdor-Kleingrothaus et al. 
reported the first evidence for the neutrinoless double beta decay [54]. However, their result was strongly 
criticized by some authors [55] . Future experiments will have sufficiently high sensitivity to clarify the 
present debates [56], either to confirm or to disprove the alleged result in Ref. [54]. 
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specific regions. To see this point more clearly, we use p and a to denote the Majorana 
CP- violating phases of V and take the convention [36] 



arg(Ki)=P, arg(y e2 )=<7, arg(K 3 ) = . (2.32) 

Of course, p and a have nothing to do with CP and T violation in normal neutrino 
oscillations. In view of (2.30) and (2.32), we find that (m) ee = requires 



mi 


\V e2 


2 


sin 2a 


m 2 


\v el 
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sin 2p 


m 2 


\Ve3 


2 


sin2p 


m 3 


\v e2 


2 


sin2(cr — p) 



(2.33) 

As rrii > holds, part of the (p, <r) parameter space must be excluded. To pin down the 
whole ranges of p and a allowed by current neutrino oscillation data under the (m) ee = 
condition, we utilize (2.18) and (2.33) to calculate the ratio of solar and atmospheric 
neutrino mass-squared differences: 



Am s 2 un \V f 
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Ki| sin 2 2p- \ V e2 \ A sm 2 2a 



Am 2 tm \V el \ 



\/ e2 | 4 sin 2 2(a-p) - |\/ e3 | 4 sin 2 2p 



(2.34) 



Then Am 2 un /Am 2 tm ~ 10~ 2 imposes quite strong constraints on the (p, a) parameter 
space, as numerically illustrated in Fig. 2 of Ref. [58]. One can in turn arrive at the 
allowed ranges of mi/m 2 and m 2 /m 3 with the help of (2.33). 

We certainly hope that (m) ee is non-zero and measurable. If the parameters of 
neutrino oscillations are well determined, a measurement of (m) ee will allow us to extract 
valuable information about the Majorana phases of CP violation and the absolute scale 
of neutrino masses. Taking account of (2.19) and (2.32), one may rewrite (2.30) as 



(m)ee = 



ym|+pAm 2 tm + gAm s 2 un \V el \ 2 e 2ip 

(2.35) 



+ J mi +pAm 2 tm \V e2 \ z e Zl ° + m 3 \V t 



e3| 



where \V ei \ can be expressed in terms of the mixing factors of solar and CHOOZ reactor 
neutrino oscillations (see section 3.1.1 for details): 

\V el \ 2 = i ^cos 2 6> chz + \J cos 4 9 chz - sin 2 26» sun ^ , 

|V e2 | 2 = - (cos 2 9 chz - \J cos 4 9 chz - sin 2 29 smi ^j , 
\V e3 \ 2 = sin 2 9 chz . (2.36) 
We see that (m) ee consists of three unknown parameters: m 3 , p and a, which are unable 



to be determined in any neutrino oscillation experiments. Once Am 2 un , Am 2 
and 9 c h z are measured to an acceptable degree of accuracy, one should be able to get 
a useful constraint on the absolute neutrino mass m 3 for arbitrary values of p and a 
from the observation of (m) ee . If the magnitude of m 3 could roughly be known from 
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Figure 2.6: Illustrative dependence of (m) ee on m 3 for the neutrino mass spectrum 
m 2 > m 3 (curves a and b) and the neutrino mass spectrum m 2 < m 3 (curves c and 
d), where we have typically taken {p, a} = {0,0} or {7r/4,7r/4} (curves a and c) and 
{p, a} = {0, 7r/4} or {7r/4, 0} (curves b and d). The region between two dashed lines 
corresponds to the values of (m) ee reported in Ref. [54]. 

some cosmological constraints, it would be likely to obtain some loose but instructive 
information on the Majorana phases p and a by confronting (2.35) with the experimental 
result of (m) ee . 

Let us illustrate the dependence of (m) ee on m 3 , p and a in a numerical way. For 
simplicity, we typically take Am^ = 5 x 10" 5 eV 2 , Am 2 tm = 3 x 10~ 3 eV 2 , q = — 1, 
sin 26 snn = 0.8 and sin 2 26 chz = 0.05 in our calculations. We consider four instructive 
possibilities for the unknown Majorana phases p and a: (1) p = a = 0; (2) p = tt/4 
and <t = 0; (3) p = and er = ir/4; and (4) p = a = 7r/4. The results for (m) ee 
as a function of m 3 are shown in Fig. 2.6. We find that it is numerically difficult to 
distinguish between possibilities (1) and (4), or between possibilities (2) and (3). The 
reason is simply that the first two terms on the right-hand side of (2.35) dominate 
over the third term, and they do not cancel each other for the chosen values of p and 
a. We also find that the changes of (m) ee are rather mild. If reasonable inputs of 
sin 2 2# sun , sin 2 2# atm and sin 2 2# c h z are taken, a careful numerical scan shows that the 
magnitude of {m) ee does not undergo any dramatic changes for arbitrary p and a. Thus 
a rough but model-independent constraint on the absolute scale of neutrino masses can 
be obtained from the observation of (m) ee . Taking account of the alleged experimental 
region 0.05 eV < (m) ee < 0.84 eV in Ref. [54], one may arrive at the conclusion that m 3 
is most likely to lie in the range 0.1 eV < m 3 < 1 eV (see Fig. 2.6). This result implies 
that mi m vri2 ~ m 3 and m\ + iri2 + m 3 pa 3m 3 pa (0.3 — 1.5) eV hold. Such a sum of 
three neutrino masses is only partly compatible with the robust WMAP upper bound 
given in (2.22). 

Because of its importance, the neutrinoless double-/? decay has attracted a lot of 
very delicate studies. For recent comprehensive analyses of the dependence of (m) ee on 
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neutrino masses and Majorana phases, we refer the reader to Refs. [43, 55, 59] and 
references therein. 



2.3 Lepton Mass Matrices 

At low energy scales, the phenomenology of lepton masses and flavor mixing can be 
formulated in terms of the charged lepton mass matrix Mi and the (effective) neutrino 
mass matrix M v . They totally involve twelve physical parameters: 

• three masses of charged leptons (m e , and m T ), which have precisely been mea- 



• three masses of neutrinos (m 1; m 2 and m^), whose relative sizes have roughly been 
known from solar and atmospheric neutrino oscillations; 

• three angles of flavor mixing (0 X , 9 y and 9 Z ), whose values have been determined 
or constrained to an acceptable degree of accuracy from solar, atmospheric and 
reactor neutrino oscillations; 

• three phases of CP violation (5, p and a), which are completely unrestricted by 
current neutrino data. 

The long-baseline neutrino oscillation experiments are expected to determine or constrain 
Am^, Am§ 2 , X , 9 y , 9 Z and 5 (the Dirac-type CP- violating phase) to a high degree 
of accuracy in the future. The precision experiments for the tritium beta decay and 
the neutrinoless double-/? decay may help determine or constrain the absolute scale 
of neutrino masses. The possible regions of two Majorana-type CP-violating phases 
(p and a) could roughly be constrained by the neutrinoless double-/? decay. It seems 
hopeless, at least for the presently conceivable sets of feasible neutrino experiments, to 
fully determine (mi,m 2 ,m 3 ) and (p, a) [60]. In this case, to establish testable relations 
between the experimental quantities and the parameters of Mi and M v has to rely on 
some phenomenological hypotheses. 

The specific structures of Mi and M v depend substantially on the flavor basis to be 
chosen. There are two simple possibilities to choose the flavor basis: (a) Mi is diagonal 
and has positive eigenvalues; and (b) M v is diagonal and has positive eigenvalues. Let 
us focus on case (a) in the subsequent discussions, because it has popularly been taken 
in the phenomenological study of lepton mass matrices 7 . 

In the flavor basis where Mi is diagonal (i.e., the flavor eigenstates of charged leptons 
are identified with their mass eigenstates), the symmetric neutrino mass matrix M v can 
be diagonalized by a single unitary transformation, 



7 Case (b) has been considered in Ref. [61], for example, to derive the democratic (nearly bi-maximal) 
neutrino mixing pattern. 



sured [6]; 
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m 2 
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(2.37) 
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Table 2.1: Nine two-zero textures of the neutrino mass matrix, which are compatible 
with current neutrino oscillation data. Note that the patterns in the same category 
(e.g., A x and A 2 ) have very similar phenomenological consequences [63]. 
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The 3x3 unitary matrix V is just the lepton flavor mixing matrix, consisting of three 
mixing angles and three CP- violating phases. It is straightforward to fix the form of M v 
from (2.37), 

/mi \ 
M v — V m 2 \V T , (2.38) 
V mj 

but extra assumptions have to be made for M v in order to relate neutrino masses to flavor 
mixing parameters. The general idea is to assume that some independent matrix elements 
of M v are actually dependent upon one another, caused by an underlying (broken) flavor 
symmetry. In particular, this dependence becomes very simple and transparent, if the 
relevant matrix elements are exactly equal to zero. Three general categories of M v with 
texture zeros have been classified and discussed 8 : 



• The three- zero textures of M v . There are totally twenty possible patterns of M v 
with three independent vanishing entries, but none of them is allowed by current 
neutrino oscillation data [62]. 

• The two-zero textures of M v . There are totally fifteen possible patterns of M v with 
two independent vanishing entries. Nine two-zero patterns of M v , as illustrated in 
Table 2.1, are found to be compatible with current experimental data (but two of 
them are only marginally allowed [63]). 

• The one-zero textures of M v . There are totally six possible patterns of M v with 
one vanishing entry, as shown in Table 2.2. It has been found in Ref. [64] that 

'Because M v is symmetric, a pair of off-diagonal texture zeros of M„ have been counted as one zero. 



23 



Table 2.2: Six possible one-zero textures of the neutrino mass matrix. They may agree 
or disagree with current neutrino oscillation data, if one mass eigenvalue (mi or m 3 ) 
vanishes [64]. 
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three one- zero textures of M v with m x = and four one-zero textures of M v with 
m 3 = are compatible with the present neutrino oscillation data. 



Of course, the specific texture zeros of the neutrino mass matrix may not be preserved 
to all orders or at any energy scales in the unspecified interactions from which neutrino 
masses are generated. Whether mi = or m 3 = is stable against radiative corrections 
should also be taken into account in a concrete theoretical model. 

In general, both the charged lepton mass matrix Mi and the (effective) neutrino mass 
matrix M v are not diagonal. A typical ansatz of this nature is the Fritzsch pattern [65] 
of lepton mass matrices, 



(2.39) 



r n in Ref. [66] that such 
a simple ansatz of lepton mass matrices can naturally predict a normal hierarchy of 
neutrino masses and a bi-large pattern of lepton flavor mixing. A generalized version of 
(2.39) with four texture zeros, 
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(2.40) 



has carefully been analyzed in Ref. [67]. The motivation to consider the four-zero texture 
of lepton mass matrices is quite clear. It is well known that the four-zero texture of 
quark mass matrices is more successful than the six-zero texture of quark mass matrices 
to interpret the strong hierarchy of quark masses and the smallness of flavor mixing 
angles [68]. Thus the spirit of lepton-quark similarity motivates us to conjecture that 
the lepton mass matrices might have the same texture zeros as the quark mass matrices. 
Such a phenomenological conjecture is indeed reasonable in some specific models of 
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grand unified theories [69], in which the mass matrices of leptons and quarks are related 
to each other by a new kind of flavor symmetry. That is why the four-zero texture of 
lepton mass matrices has been considered as a typical example in a number of interesting 
model-building works [70]. 

Taking account of the seesaw mechanism, one may investigate how the texture zeros 
of M D and M R are related to those of M v via (2.9). For example, it has been shown in 
Ref. [71] that some two- zero textures of M v can easily be reproduced from Md and Mr 
with several vanishing entries. Provided both Md and Mr have texture zeros in their 
(1,1), (1,3) and (3,1) positions, M„ must have the same texture zeros [12]. This seesaw 
invariance of lepton mass matrices implies that all lepton and quark mass matrices might 
have a universal texture like (2.40), resulting from a universal (broken) flavor symmetry 
at the energy scale where the seesaw mechanism works. 

There are many other phenomenological ansatze of lepton mass matrices, either at 
low energy scales or at high energy scales. For recent reviews with extensive references, 
we refer the reader to Ref. [12], Refs. [18]-[21], and Ref. [69]. 

3 Lepton Flavor Mixing 

3.1 The Flavor Mixing Matrix in Vacuum 

The phenomenon of lepton flavor mixing arises from the mismatch between the diagonal- 
ization of the charged lepton mass matrix M\ and that of the neutrino mass matrix M v 
in an arbitrary flavor basis. Without loss of any generality, one can choose to identify 
the flavor eigenstates of charged leptons with their mass eigenstates. In this specific 
basis with M\ being diagonal, M v is in general not diagonal. If massive neutrinos are 
Dirac particles, it is always possible to diagonalize the arbitrary mass matrix M v by a 
bi-unitary transformation: 



If massive neutrinos are Majorana particles, one may diagonalize the symmetric mass 
matrix M v by a single unitary transformation: 

/m, \ 



In (3.1) and (3.2), m i denote the (positive) mass eigenvalues of three active neutrinos. 
In terms of the mass eigenstates of charged leptons and neutrinos, the Lagrangian of 
charged-current weak interactions in (2.1) can be rewritten as 




(3.1) 




(3.2) 



9 



K^^vVj^J w; + h.c. , 




(3.3) 
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where the 3x3 unitary matrix V links the neutrino mass eigenstates v 2 ,v ? ) to the 
neutrino flavor eigenstates (u e , u T ): 





fv el 


v e2 


v e3 \ 






= 










VKi 


v t2 


v T3 J 





Obviously V is just the flavor mixing matrix of charged leptons and active neutrinos. 
The analogue of V in the quark sector is well known as the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [10]. In some literature, the lepton mixing matrix V has been referred 
to as the Maki-Nakagawa-Sakata (MNS) matrix [27] or the Pontecorvo-MNS (PMNS) 
matrix. 

Note that massive neutrinos are commonly believed to be Majorana particles, al- 
though there has not been any reliable experimental evidence for that. For simplicity, 
we use V to denote the lepton flavor mixing matrix in the Majorana case. Then V can 
always be parametrized as a product of the Dirac-type flavor mixing matrix U and a 
diagonal phase matrix P; i.e., V — UP [72]. The nontrivial CP- violating phases in P 
characterize the Majorana nature of V. It is straightforward to show that the rephasing- 
invariant quantities V^V^V^V^ and U ai U pjU^U^ are identical for arbitrary indices 
(cx,(3) over (e, //, r) and over (1,2,3). This result implies that V and U have the 
same physical effects in normal neutrino-neutrino and antineutrino-antineutrino oscilla- 
tions, whose probabilities depend only upon the rephasing-invariant quantities mentioned 
above. One then arrives at the conclusion that it is impossible to determine the nature of 
massive neutrinos by studying the phenomena of neutrino oscillations. Only the exper- 
iments which probe transitions between left-handed and right-handed neutrino states, 
like the neutrinoless double beta decay, could tell whether massive neutrinos are of the 
Majorana type or not. 



3.1.1 Experimental Constraints 

Current experimental constraints on the lepton flavor mixing matrix V come mainly 
from solar, atmospheric and CHOOZ reactor neutrino oscillations. These three types of 
neutrino oscillations all belong to the "disappearance" experiments, in which the survival 
probability of a flavor eigenstate v a is given as 

P(u a ^u a ) = 1 - 4X;(|K;| 2 |K/sin 2 i^) , (3.5) 

i<j 

where = \ .lHS.vri 2 - i L j E with Ara^ = m| — mf, L stands for the baseline length 
(in unit of km), and E is the neutrino beam energy (in unit of GeV). Because of CPT 
invariance, the survival probabilities of neutrinos and antineutrinos in vacuum are equal 
to each other: P(v a — > V a ) = P{y a — > u a ). This equality may in general be violated, if 
neutrinos and antineutrinos propagate in matter [73]. 

As mentioned in section 2, the solar {y e — > v e ) and atmospheric iy^ — > v^) neutrino 
oscillation data indicate 

\Am 2 21 \ = Am 2 SUQ < Am^ tm = |Am^ 2 | » \Am 2 31 \ . (3.6) 
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In addition, the CHOOZ {y e — > u e ) neutrino oscillation data are obtained at the scale 
Am 2 hz S3 |Am| 2 |. To a good degree of accuracy, (3.5) is simplified to 

/ Am 2 T \ 

P{v e - u e ) S3 1 -sin 2 2£ sun sin 2 f 1.27 ^ un j (3.7) 

with 

sin 2 2fl sun = A\V el \ 2 \V e2 \ 2 (3.8) 

for solar neutrino oscillations; 

(Am 2 T \ 
1.27 | m j (3.9) 

with 

sin 2 20 atm = 4|V^ 3 | 2 (l - |^ 3 | 2 ) (3.10) 
for atmospheric neutrino oscillations; and 

/ Am 2 T \ 

P(V e ^V e ) S3 1 - sin 2 29 chz sin 2 f 1.27 ^ j (3.11) 

with 

sin 2 2£ chz = 4|K 3 | 2 (l - IK3I 2 ) (3.12) 

for the CHOOZ neutrino oscillation experiment. In view of the present experimental 
data from SK [2], SNO [3], KamLAND [4], K2K [5] and CHOOZ [7], we have 



0.25 < sin 2 sun < 0.4 or 30.0° < sun < 39.2° , 



2 , 

0.92 < sin 2 2# atm < 1.0 or 36.8° < atm < 53.2°, 

< sin 2 2# chz < 0.1 or 0° < 9 chz < 9.2°, (3.13) 

at the 90% confidence level. More precise data on 9 sun , atm and 9 chz will be available in 
the near future. 

Taking account of the unitarity of V, one may reversely express \V e i\, |Vy, | V^3 1 , 
|V^ 3 | and \V t3 \ in terms of sun , 6> atm and 9 chz [36]: 



\V e i\ = \/ cos 2 6 chz + V cos 4 6 chz - sin 2 2# sun , 



1 
1 

V2 



\V e2 \ = -^= V cos2 ^chz - V cos 4 6 chz - sin 2 26» sun , 



\V e3 \ = sin0 ch z , 
| | = sin0 atm , 

|Ks| = \l cos 2 # chz - sin 2 9 atm . (3.14) 

The other four matrix elements of V (i.e., |V^i|, |V^ 2 |, |V^i| and | V^- 2 1 ) are entirely un- 
restricted, however, unless one of them or the rephasing invariant of CP violation of V 
(defined in section 4.1) is measured. A realistic way to get rough but useful constraints 
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on those four unknown elements is to allow the Dirac phase of CP violation in V to 
vary between and 7r [74], such that one can find out the maximal and minimal mag- 
nitudes of each matrix element. To see this point more clearly, we adopt a simplified 
parametrization of V in which the Majorana phases of CP violation are omitted [75] : 



V | C x SyS Z 



S X SyS Z ~\~ C X CyC 



-iS 



~$xCy$z 



x^y^ 



C X SyC 



SyC Z 



CyC Z 



(3.15) 



where s x = sin 6^, c x = cos9 x , and so on. The Dirac phase 5 affects the magnitudes 
of V^i, V^ 2 , V T i and V r2 , while the Majorana phases of CP violation do not have such 
effects (see section 3.1.2 for more discussions). Note that three mixing angles (9 X , 9 yi 9 Z ), 
which are all arranged to lie in the first quadrant, can be written as 
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(3.16) 



It is then straightforward to obtain 
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(3.17) 



Varying the Dirac phase 5 from to n, we are led to the most generous ranges of |V^i|, 
|V^ 2 |, |Ki| and |K 2 | [76]: 
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(3.18) 



Note that the lower and upper bounds of each matrix element turn to coincide with each 
other in the limit \V e3 \ — > 0. Because of the smallness of \V e3 \, the ranges obtained in 
(3.18) should be quite restrictive. Hence it makes sense to recast the lepton flavor mixing 
matrix even in the absence of any experimental information on CP violation. 

Using the values of 9 sun , 6> atm and 9 c h z given in (3.13), we calculate \V e i\, \V e2 \, \V e3 \, 
\V^ 3 \ and \V t3 \ from (3.14). Then the allowed ranges of |V^i|, \V ll2 \ ) \V T i\ and \V r2 \ can 
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be found with the help of (3.18). Our numerical results are summarized as [76] 

/0.70-0.87 0.50 -0.69 < 0.16 
\V\ = 0.20 - 0.61 0.34- 0.73 0.60 - 0.80 
\0.21 - 0.63 0.36 - 0.74 0.58 - 0.80 

This represents our present knowledge on the lepton flavor mixing matrix. 



(3.19) 



3.1.2 Standard Parametrization 



Flavor mixing among n different lepton families is in general described by a n x n unitary 
matrix V, whose number of independent parameters relies on the style of neutrinos. If 
neutrinos are Dirac particles, V can be parametrized in terms of n(n — l)/2 rotation 
angles and (n — l)(n — 2)/2 phase angles. If neutrinos are Majorana particles, however, 
a full parametrization of V requires n(n — l)/2 rotation angles and the same number of 
phase angles 9 . The flavor mixing of charged leptons and Dirac neutrinos is completely 
analogous to that of quarks, for which a number of different parametrizations have been 
proposed and classified in the literature [78]. In this subsection we classify all possible 
parametrizations for the flavor mixing of charged leptons and Majorana neutrinos with 
n — 3. Regardless of the phase- assignment freedom, we find that there are nine struc- 
turally different parametrizations for the 3x3 lepton flavor mixing matrix V. Although 
all nine representations are mathematically equivalent, one of them is likely to make the 
underlying physics of lepton mass generation (and CP violation) more transparent, or is 
more useful in the analyses of neutrino oscillation data. 

The 3x3 lepton flavor mixing matrix V defined in (3.4) can be expressed as a product 
of three unitary matrices 0\, 2 and 3 , which correspond to simple rotations in the 
complex (1,2), (2,3) and (3,1) planes: 
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where Sj = sin^j and q = cos^j (for i — 1,2, 3). Obviously OiO\ = 0\Oi = 1 holds, and 
any two rotation matrices do not commute with each other. We find twelve different 
ways to arrange the product of 0\, 2 and 3 , which can cover the whole 3x3 space 
and provide a full description of V. Explicitly, six of the twelve different combinations 
of Oi belong to the type [79] 

V = Oije^P^h ) ® 0,(9,, a 3 , (3j, 7j ) ® O t (6[, a[, (3[, 7 (3.21) 

9 No matter whether neutrinos are Dirac or Majorana particles, the n x n flavor mixing matrix has 
(n — l) 2 (n — 2) 2 /4 rephasing invariants of CP violation [77]. 
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with % ^ j, where the complex rotation matrix 0« occurs twice; and the other six belong 
to the type 

V = O0 t , a^Pi, 7l ) ® Oj(6j, a^Pj^j) ® O k (9 k , a k , [3 k , lk ) (3.22) 

with i 7^ j 7^ k, in which the rotations take place in three different complex planes. Note 
that the products OiOjOi and 0{O k Oi (for % ^ k) in (3.21) are correlated with each 
other, if the relevant phase parameters are switched off. Hence only nine of the twelve 
parametrizations, three from (3.21) and six from (3.22), are structurally different. 

In each of the nine distinct parametrizations for V, there apparently exist nine phase 
parameters. Six of them or their combinations can be absorbed by redefining the arbi- 
trary phases of charged lepton fields and a common phase of neutrino fields. If neutrinos 
are Dirac particles, one can also redefine the arbitrary phases of Dirac neutrino fields, 
so as to reduce the number of the remaining phase parameters from three to one. In 
this case V contains a single CP-violating phase of the Dirac nature. If neutrinos are 
Majorana particles, however, there is no freedom to rearrange the relative phases of Ma- 
jorana neutrino fields. Hence three nontrivial phase parameters are present in the 3x3 
Majorana-like flavor mixing matrix V. Two of the three CP-violating phases in V can 
always be factored out as the "pure" Majorana phases through a proper phase assign- 
ment of the charged lepton fields, and the remaining one is identified as the "Dirac-like" 
phase. Both CP- and T-violating effects in normal neutrino oscillations depend only 
upon the Dirac-like phase. 

Of course, different parametrizations of V are mathematically equivalent, and adopt- 
ing any of them does not point to physical significance. It is quitely likely, however, 
that one particular parametrization is more useful and transparent than the others in 
the analyses of data from various neutrino experiments and (or) towards a deeper un- 
derstanding of the underlying dynamics responsible for lepton mass generation and CP 
violation. We find that such an interesting parametrization does exist [79]: 

^ Cjj; S X C g t 
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C x SyS z S x CyC ' S X SyS Z ~\~ C X CyC ' ^y^Z 








0\ 
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1/ 



(3.23) 



with s x = sin 61,,, c y = cos 6^, and so on. Without loss of generality, the three mixing 
angles (9 X , 9 y , 6 Z ) can all be arranged to lie in the first quadrant. Arbitrary values between 
and 2tc are allowed for the Dirac CP-violating phase 5 and the Majorana CP-violating 
phases p and a. 

The parametrization in (3.23) can be regarded as a straightforward generalization 
of the parametrization in (3.15) to include the Majorana phases of CP violation. Thus 
the interesting relations obtained in (3.16) remain valid. A remarkable merit of this 
parametrization is that its three mixing angles (9 X , 9 y , 9 Z ) are directly related to the 
mixing angles of solar, atmospheric and CHOOZ reactor neutrino oscillations: 

9 X ~ 9 sml , 9 y ~ ^atm , 9 Z ~ ^ c h z , (3.24) 

which can easily be derived from (3.8), (3.10), (3.12) and (3.23). Other parametrizations 
of the lepton flavor mixing matrix would not be able to result in such simple relations 
between the fundamental mixing quantities and the relevant experimental observables. 



30 



The CP- and T-violating effects in normal neutrino oscillations are measured by the 
well-known Jarlskog parameter J [80] (defined in section 4.1), which is proportional to 
sin 5 in the following way: 

J ps sin 9 smi cos 9 smi sin atm cos 6» atm sin 9 chz cos 2 9 chz sin 5 . (3.25) 

To see this point more clearly, we write out the transition probabilities of different 
neutrino flavor eigenstates in vacuum [12]: 

P(v a ^v p ) = -^[^(v^V^V^sin 2 ^,] - 8 J C 8 ^ ^"0 > ( 3 - 26 ) 

i<j i<j 

where (cx,f3) run over (e, //), (//, r) or r, e); i and j run over 1, 2, 3; and Fji and Am^ 
have been defined below (3.5). In the assumption of CPT invariance, the transition 
probabilities P{y$ — > v a ) or P(V a — > 77^) can directly be read off from (3.26) through 
the replacement J ==>- —J (i.e., V ==>- V*). Therefore, the CP-violating asymmetry 
between P{y a — > Vp) and P(V a — > I7g) amounts to the T-violating asymmetry between 
P{y a -> i/^) and P(z/ /3 -> i/ a ): 

AP = P(u a ^u p )-P(u a ^u p ) 

= — 16Jsin P 2 i sin P 31 sinP 32 . (3.27) 

The determination of J from AP will allow us to extract the CP-violating phase 5, 
provided all three mixing angles (9 x ,9 y ,9 z ) have been measured elsewhere. In practical 
long-baseline neutrino oscillation experiments, however, all these measurable quantities 
may be contaminated by the terrestrial matter effects. Hence the fundamental parame- 
ters of lepton flavor mixing need be disentangled from the matter-corrected ones. Some 
discussions about the matter effects on lepton flavor mixing and CP violation will be 
given in sections 3.2 and 4.1. 

Regardless of the Majorana phases p and a, which have nothing to do with normal 
neutrino oscillations, we have located the Dirac phase 5 in such a way that the matrix 
elements in the first row and the third column of V are real. As a consequence, the CP- 
violating phase 5 does not appear in the effective mass term of the neutrinoless double 
beta decay. Indeed, the latter reads: 

(m) ee = mil/ 2 + m 2 V e 2 2 + m 3 V^ 



where 



= A/a+ bcos2(p — cr) + ccos2p + dcos2a , (3.28) 



2 4 4, 244, 24 

a = m 1 c x c z + m 2 s x c z + m 3 s z 
b = 2mim 2 s 2 c 2 c^ , 

o 2 2 2 

c = 2mim 3 c x s z c z , 



d = 2m 2 m 3 s x s;c 2 z . (3.29) 

It becomes obvious that (m) ee is independent of the Dirac phase 5. On the other hand, 
the CP- or T-violating asymmetry AP in normal neutrino oscillations is independent of 
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the Majorana phases p and o. Thus two different types of CP-violating phases can (in 
principle) be studied in two different types of experiments. 

A long-standing and important question is whether the two Majorana phases p and a 
can separately be determined by measuring other possible lepton-number-nonconserving 
processes, in addition to the neutrinoless double beta decay. While the answer to this 
question is affirmative in principle, it seems to be negative in practice. The key problem 
is that those lepton-number-violating processes, in which the Majorana phases can show 
up, are suppressed in magnitude by an extremely small factor compared to normal weak 
interactions [79, 81]. Therefore it is extremely difficult, even impossible, to measure or 
constrain p and a in any experiment other than the one associated with the neutrinoless 
double beta decay. 

On the theoretical side, how to predict or calculate those flavor mixing angles and CP- 
violating phases on a solid dynamical ground is also a big challenge. Phenomenologically, 
the parametrization in (3.23) is expected to be very useful and convenient, and might even 
be able to provide some insight into the underlying physics of lepton mass generation. We 
therefore recommend it to experimentalists and theorists as the standard parametrization 
of the 3x3 lepton flavor mixing matrix. 

3.1.3 Constant Mixing Patterns 

Combining (3.13) and (3.24), one can immediately observe the essential feature of lepton 
flavor mixing: two mixing angles are large and comparable in magnitude (6 X ~ y ~ 1), 
while the third one is very small (9 Z <C 1). Such a "nearly bi-maximal" flavor mixing 
pattern has attracted a lot of attention in model building. In particular, a number of 
constant mixing patterns, which are more or less close to the experimentally-allowed 
lepton flavor mixing matrix in (3.19), have been proposed. Possible flavor symmetries 
and their spontaneous or explicit breaking mechanisms hidden in those constant patterns 
might finally help us pin down the dynamics responsible for neutrino mass generation 
and lepton flavor mixing. It is therefore worthwhile to have an overview of a few typical 
constant lepton mixing patterns existing in the literature 10 . 

The first example is the so-called "democratic" lepton flavor mixing pattern, 
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(3.30) 



which was originally obtained by Fritzsch and Xing [61] from the breaking of flavor 
democracy of the charged lepton mass matrix in the flavor basis where the neutrino 
mass matrix is diagonal. The predictions of (3.30) for the mixing factors of solar and 

10 We only concentrate on those nearly bi-maximal flavor mixing patterns, in which 6 X <~ 9 y >> 6 Z 
holds. Hence the interesting Cabibbo ansatz [82] , which has \V a i 1 = 1/ V3 (for a = e,fi,T and i = 1, 2, 3) 
or 6 X = By = 45° and 9 Z w 35.3°, will not be taken into account. Such a "tri-maximal" mixing pattern 
has completely been ruled out by current experimental data. 
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Table 3.1: The mixing angles (0 x ,9 y ,9 z ) in four constant lepton flavor mixing patterns. 



Mixing Angle 


Pattern (3.30) 


Pattern (3.31) 


Pattern (3.32) 


Pattern (3.33) 


9 X 


45° 


45° 


35.3° 


30° 


9y 


54.7° 


45° 


45° 


45° 


e z 


0° 


0° 


0° 


0° 



atmospheric neutrino oscillations, sin 2 29 sun = 1 and sin 2 2# atm = 8/9, are not favored 
by today's experimental data. Hence slight modifications of (3.30) have to be made. In 
Ref. [61] and Ref. [83], some proper perturbative terms have been introduced to break 
the S(3)l x S(3)r symmetry of the charged lepton mass matrix and the S(3) symmetry 
of the neutrino mass matrix. The resultant flavor mixing matrix, whose leading term 
remains to take the form of (3.30), is able to fit current solar, atmospheric and CHOOZ 
reactor neutrino oscillation data very well. 

The second example is the so-called "bi-maximal" neutrino mixing pattern, 



V = 



proposed by Barger et al [84] and by Vissani [85]. It predicts sin 2 29 s 



( V2 
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V2 
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1 
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V \ 
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2 


V2 
2 / 



(3.31) 



sm J 26» atm = 1 

for solar and atmospheric neutrino oscillations. As the result sin 2 20 sun = 1 is disfavored 
by the LMA solution to the solar neutrino problem, slight modifications of (3.31) are also 
required. A number of possibilities to modify the bi-maximal mixing pattern into nearly 
bi-maximal mixing patterns have been discussed in the literature [86]. The essential 
point is to introduce a perturbative term to (3.31), which may naturally be related to 
the mass ratio (s) of charged leptons or to the Cabibbo angle of quark flavor mixing, such 
that sin 2 29 sun deviates from unity to a proper extent. 

Another interesting neutrino mixing pattern is 



(3.32) 



It was first conjectured by Wolfenstein long time ago [87] 11 . Possible (broken) flavor 
symmetries associated with (3.32) have recently been discussed by several authors [88]. 
Two straightforward consequences of (3.32) on neutrino oscillations are sin 2 26 l sun = 8/9 

11 The first and second columns of V were interchanged in the original paper of Wolfenstein [87] . Such 
an interchange is allowed by current neutrino oscillation data, because \ V e i\ and \ V e 2\ cannot separately 
be determined from the measurement of sin 2 29 SUD 
a common convention to set > \V e 2\- 
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4|V;i| 2 |V;2r- In (3.14) and (3.32), we have taken 
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and sin 2 2# atm = 1, favored by the present solar and atmospheric neutrino oscillation 
data. Of course, nonvanishing V e 3 (or 9 Z ) and CP-violating phases can be introduced 
into (3.32), if proper perturbations to the charged lepton and neutrino mass matrices 
are taken into account. 

The last example for constant lepton flavor mixing is given as 
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V = 




V2 
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V2 
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%/2 
2 / 



which has recently been discussed by Giunti [89], Peccei [90] and Xing [91]. It predicts 
sin 2 29 sun = 3/4 and sin 2 2# atm = 1, which are in good agreement with the present solar 
and atmospheric neutrino oscillation data. Again slight modifications of (3.33) can be 
done [91], in order to accommodate nonvanishing V e 3 and leptonic CP violation. 

We summarize the values of three flavor mixing angles for each of the four constant 
patterns in Table 3.1. The similarity and difference between any two patterns can easily 
be seen. It is worth remarking that a specific constant mixing pattern should be regarded 
as the leading-order approximation of the "true" lepton flavor mixing matrix, whose 
mixing angles depend in general on both the ratios of charged lepton masses and those 
of neutrino masses. Naively, one may make the following speculation: 

• Large values of 9 X and 9 y could arise from a weak hierarchy or a near degeneracy 
of the neutrino mass spectrum, as the strong hierarchy of charged lepton masses 
implies that m e /m^ and m^/m T are unlikely to contribute to 9 X and 9 y dominantly. 

• Special values of 9 X and 9 y might stem from a discrete flavor symmetry of the 
charged lepton mass matrix or the neutrino mass matrix. Then the contributions 
of lepton mass ratios to flavor mixing angles, due to flavor symmetry breaking, are 
only at the perturbative level. 

• Vanishing or small 9 Z could be a natural consequence of the explicit textures of 
lepton mass matrices. It might also be related to the flavor symmetry which gives 
rise to sizable 9 X and 9 y . 

• Small perturbative effects on a constant flavor mixing pattern can also result from 
the renormalization-group equations of leptons and quarks, e.g., from the high 
energy scales to the low energy scales or vice versa [92]. 

To be more specific, we refer the reader to Ref. [12] and Ref. [21], from which a lot of 
theoretical models and phenomenological ansatze can be found for the interpretation of 
possible spectra of lepton masses and possible patterns of lepton flavor mixing. 

3.2 The Flavor Mixing Matrix in Matter 

The effective Hamiltonian responsible for the propagation of neutrinos in vacuum can be 
written as 7i eff = Q u / (2E) in the basis where the flavor eigenstates of charged leptons 
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are identified with their mass eigensates. In Ti e s : E m i denotes the neutrino beam 
energy, and Q u is given by 

(m\ \ 
0, = V m 2 2 \V ] , (3.34) 
V m 2 J 

in which V is the 3x3 lepton flavor mixing matrix. Taking account of (3.1) or (3.2), we 
obtain Vl u = M v Ml [93], no matter whether massive neutrinos are Dirac or Majorana 
particles. A simple relation can therefore be established between the effective Hamil- 
tonian 7i c ^ and the neutrino mass matrix M u . The form of 7i C R has to be modified, 
however, in order to describe the propagation of neutrinos in matter. 

If neutrinos travel through a normal material medium (e.g., the earth), which consists 
of electrons but of no muons or taus, they encounter both charged- and neutral-current 
interactions with electrons. In this case, the effective Hamiltonian responsible for the 
propagation of neutrinos takes the form = (2E) , where [33] 

(A 0\ 

n„ = n„ + ooo. (3.35) 
V o oo; 

The apparent deviation of Q u from Q u arises from the matter-induced effect; namely, A = 
2\^2 GfN £ E describes the charged-current contribution to the u e e~ forward scattering, 
where N e is the background density of electrons and E stands for the neutrino beam 
energy. The neutral-current contributions, which are universal for u e , and v T neutrinos 
and can only result in an overall unobservable phase, have been omitted. 

Let us concentrate on the long-baseline neutrino oscillation experiments, which are 
designed to measure the parameters of lepton flavor mixing and CP violation. We shall 
assume a constant earth density profile (i.e., N e = constant) in section 3.2.1 and section 
3.2.3. Such an assumption is rather plausible and close to reality, only if the neutrino 
beam does not go through the earth's core [94] . 

In the chosen flavor basis where Mi is diagonal, VL V can be diagonalized by a unitary 
transformation: 

/m? \ 

V ] , (3.36) 

where rhi are the effective neutrino masses in matter. The unitary matrix V is just 
the matter-modified lepton flavor mixing matrix. Denoting the effective neutrino mass 
matrix in matter as M u , one may get Cl u = \ 1 1 similar to Q u = \ /,.. Phenomeno- 
logically, it is important to find the analytical relations between and m 8 as well as 
the exact relation between V and V. 
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3.2.1 Generic Analytical Formulas 



Combining (3.34), (3.35) and (3.36), we calculate in terms of and the matrix 
elements of V. The results are 



35 
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ml 



= m, H — x 

1 3 3 



m\ + -x + -z\Jx 2 



3y 



V 3 (1 " * 2 ) 



3y , 



(3.37) 



where x, y and z are given by [95] 

x = Am^ + Am^ + A , 
y = Am 2 21 Am 2 31 + A [Am^ (l - |K 



'e2 



(1 - IV.sl 2 ) 

2" 



COS 



1 2x 3 - 9xy + 27AAm 2 Ami, I V el 
- arccos 



2 (a; 2 - 3yf 2 



(3.38) 



with Am\ x 



m n 



m\ and Am^ 



m 2 . Of course, m 2 = m\ can be reproduced 



from (3.37) when A = is taken. Only the mass-squared differences Am^ = m\ — m\ 
m 2 , which depend on Am^, Aro^, A and |V e j| 2 (for % = 1,2,3), are 



and Amj! 



~ 2 



m 



relevant to the practical neutrino oscillations in matter. 

It is easy to prove that the sum of m 2 is related to that of m 2 through 

3 3 

E™ 2 = E™ 2 + A - 

i=i i=i 



(3.39) 



Also m 2 , |Ki| 2 , A- and mf are correlated with one another through another instructive 
equation: 



A (m 2 - mfj (m 2 - mfj \V ek \ 2 = J[ (m 2 n - mfj , 



(3.40) 



n=l 



where (i,j,k) run over (1,2,3) with i ^ j ^ k. Assuming jAm^l ^> jAm^] and taking 
k — 3, one can use (3.40) to derive the approximate analytical result for the correlation 
between m 2 and mf obtained in Ref . [96] . 

The analytical relationship between the matrix elements of V and those of V can 
also be derived from (3.34), (3.35) and (3.36). After some lengthy but straightforward 
calculations [97], we arrive at 



Ni 



A 



Vai ^ Vai -r^ Ki 



D 



A 



(m 2 - mfj V e * k V ak + (m 2 - m 2 k ) V*V aj 



(3.41) 



where a runs over (e, /i, r) and k) run over (1, 2, 3) with i ^ j ^ k, and 

Ni = (m 2 -m 2 )(m 2 -m 2 )-A[(m 2 -m 2 )|K fe | 2 + (m 2 -m 2 )|K,| 2 



D 2 = N 2 + A 2 \V ei 



(in 2 - mfj \V ek \ 2 + (m 2 - mfj \V ej \ 



(3.42) 



Obviously, A = leads to V ai = V ai . This exact and compact formula shows clearly 
how the flavor mixing matrix in vacuum is corrected by the matter effects. Instructive 
analytical approximations can be made for (3.41), once the spectrum of neutrino masses 
is experimentally known or theoretically assumed. 

The results obtained above are valid for neutrinos interacting with matter. As for 
antineutrinos, the corresponding formulas for the effective neutrino masses and flavor 
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mixing matrix elements in matter can straightforwardly be obtained from (3.37)-(3.42) 
through the replacements V =>- V* and A =>- —A 

For illustration, we carry out a numerical analysis of the terrestrial matter effects on 
lepton flavor mixing and CP violation for a low-energy (100 MeV < E < 1 GeV) and 
medium-baseline (100 km < L < 400 km) neutrino oscillation experiment [93]. In view 
of current data, we typically take Am^i = 5 x 10" 5 eV 2 and Am^ = 3 x 10~ 3 eV 2 as 
well as 9 X p=* 35°, 9 y ps 40°, 9 Z ps 5° and 5 ps ±90° in the standard parametrization of 
V. In addition, the dependence of the terrestrial matter effect on the neutrino beam 
energy is given as A = 2.28 x 10" 4 eV 2 £/[GeV] [94]. With the help of (3.37) and (3.41), 
we are then able to compute the ratios Am-j/Am-j (for i — 2,3) and |V^,j|/|V^j| (for 
a = e,fi,r and i = 1,2,3) as functions of E. The relevant numerical results are shown 
respectively in Fig. 3.1 and Fig. 3.2. We see that the earth-induced matter effect on 
Awijj is significant, but that on Am| x is negligibly small in the chosen range of E. The 
matter effects on |V^3| and | j are negligible in the low-energy neutrino experiment. 
The smallest matrix element |V^| is weakly sensitive to the matter effect. In contrast, 
the other six matrix elements of V are significantly modified by the terrestrial matter 
effects. 



3.2.2 Commutators and Sum Rules 

If the charged lepton mass matrix Mi and the neutrino mass matrix M v could simulta- 
neously be diagonalized in a given flavor basis, there would be no lepton flavor mixing. 
A very instructive measure of the mismatch between the diagonalization of M\ and that 
of M v in vacuum should be the commutators of lepton mass matrices [93], defined as 





M v Ml 


= iX , 


Mi Mi 


M\M V 


EE iX' 


'M\Mi 


M v Ml 


EE iY, 


M\Mi 


M\M V 


EE iY' . 



(3.43) 
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Figure 3.2: Ratios IV^I/IV^] (for a = e,fj,,r and i = 1,2,3) changing with the beam 
energy E (in unit of GeV) for neutrinos {y) and antineutrinos (u), in which Am^ = 
5 x 10" 5 eV 2 , Am 2 31 = 3 x 10" 3 eV 2 , 9 X w 35°, 9 y w 40°, ^ 2 w 5° and 5 w ±90° have 
typically been input. 



As for neutrinos propagating in matter, the corresponding commutators of lepton mass 
matrices are defined as 



mm) 



iivii , M V M\ 

[m x m\ , .\ /;;.!/,, 

\M\Mx , M V M\ 

\m\m x , .\ /;;.!/,. 



= ix', 
= iy, 
= if. 



(3.44) 



It is obvious that these commutators must be traceless Hermitian matrices. Without loss 
of generality, one may again choose to identify the flavor eigenstates of charged leptons 
with their mass eigenstates. In this specific basis, where M t takes the diagonal form 
D t = Diag{m e , m M , m T }, we obtain X = Y, X' = Y' and X = Y , X 1 = Y' . 
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Taking account of (3.1), (3.2), (3.34) and (3.36), we obtain 

M U M\ = YDlY = o, , 

M v Ml = vblv^ = n u 



(3.45) 



in the chosen flavor basis, where D v = Diagjmi, m 2 , m 3 } and D u = Diagjmx, m 2 , m 3 }. 
In contrast, there is no direct relationship between M\M U (or MlM v ) and fJj, (or Cl v ). 
Hence we are only interested in the commutators X and X, whose explicit expressions 
read as follows: 
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L -*fLT ZJ T(l 








(3.46) 



where A a/3 = — m| for a ^ (5 running over (e, //, r), and 



= e K 2 k^*) , 
t=i 

Zap = E (^K^*) • 



(3.47) 



i=i 



One can see that A^q, = — A a @, Z/3 a = and Zp a = ZJ*^ hold. To find out how Z a/ 3 
is related to Z a/ 3, we need to establish the relation between X and X. With the help of 
(3.35), it is easy to prove that X and X are identical to each other: 



X = 



= % 



a, , Dt 



= x . 



(3.48) 



This interesting result, whose validity is independent of the specific flavor basis chosen 
above, implies that the commutator of lepton mass matrices defined in vacuum is in- 
variant under matter effects. As a straightforward consequence of X = X, we arrive at 
Z a/ 3 = Z a/ 3 from (3.46). Then a set of concise sum rules of neutrino masses emerge [98]: 



E (rfvjfa) = E y<v m v;) . 



(3.49) 



Although we have derived these sum rules in the three- neutrino mixing scheme, they may 
simply be generalized to hold for an arbitrary number of neutrino families. In addition, 
we emphasize that the useful results obtained above are independent of the profile of 
matter density (i.e., they are valid everywhere in the sun or in the earth). 

It should be noted that Z a p and Z af3 are sensitive to a redefinition of the phases 
of charged lepton fields. The simplest rephasing-invariant equality is \Z a p\ = \Z a p\, 
where (a,/3) may run over (e, //), (//, r) or (r, e). Another useful rephasing-invariant 
relationship is 

Ze[iZfj, T Z Te = Z^Z^Zt-c . (3.50) 
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As one can see in section 4.1, the imaginary parts of Z e[ii Z^ T Z Te and Z e/1 Z MT Z re are related 
respectively to leptonic CP violation in vacuum and that in matter. 

It should also be noted that three additional sum rules can similarly be derived from 
(3.35). We find that Qf = ttf holds, if and only if (a, (3) ^ (e,e). Thus the following 
equalities can be obtained: 

n e jn T / - opn? = rt e jn T / - qtvii t , 

Q^Qt T -Q^ T Qt e = Q^ e ni T - {Ql e + A) . (3.51) 
By use of (3.51) and the unitarity of V, one may arrive at 

E {rnjm)V ek V; k ) = £ (m*m>V ek V; k ) , 

i^j^k i^j^k 

£ (m\m% k V; k ) = £ (mlm)V ek V T \) , 

E (rn\™%kV; k ) = E [{m 2 im )-Am 2 k )v, k V: k \ . (3.52) 

These relations, similar to those obtained previously in Ref. [99], imply that KnVg* (for 
a ^ (3) can be expressed in terms of V^Vg* , Am^ and A [100]. In particular, we are able 
to get the analytical formula for the sides of three leptonic unitarity triangles in matter 
(see section 4.2 for the detail). 

Let us remark that the results obtained above are only valid for neutrinos propa- 
gating in vacuum and in matter. As for antineutrinos, the corresponding formulas can 
straightforwardly be written out from (3.46)-(3.52) through the replacements V V* 
and A =>- -A. 



3.2.3 Effective Mixing Parameters 

Although the present non-accelerator neutrino experiments have yielded some impressive 
constraints on three lepton flavor mixing angles (9 X , 9 y , and 9 Z ), a precise determination 
of them and a measurement of the Dirac-type CP-violating phase 5 have to rely on a 
new generation of accelerator experiments with very long baselines, including the possible 
neutrino factories. In such long- or medium-baseline neutrino experiments the terrestrial 
matter effects, which may deform the oscillating behaviors of neutrinos in vacuum and 
even fake the genuine CP- violating signals, must be taken into account. It is therefore 
desirable to figure out how the standard parametrization in (3.23) will be modified by 
the matter effects, in particular in the case of a constant earth density profile. 

The analytical relationship between the elements of the effective flavor mixing matrix 
V in matter and those of the fundamental flavor mixing matrix V in vacuum has been 
given in (3.41). Adopting the parametrization of V in (3.23), we explicitly obtain 

N- A 

Vai = lj~Vai + -JT E (T ak P k i) , (3.53) 

i i k 
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where the expressions of iVj and Di have been given in (3.42), A denotes the matter 
parameter, P = Diag{e* p , e tcr , 1} is the diagonal matrix of Majorana phases, and the 
matter-associated matrix elements T Q j read as follows [101] 12 : 
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(3.54) 



V a i holds in the case of A = 0. 

The results given in (3.53) and (3.54) indicate that the diagonal Majorana phase 
matrix P on the right-hand side of V is not affected by the matter effect. As a natural 
consequence, the effective flavor mixing matrix V in matter can be parametrized in the 
same form as V: 
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(3.55) 



with s x = sm9 x and c x = cos9 x , and so on. It should be noted that the matrix elements 
V^3 and V r3 in (3.53) are complex and dependent on the Dirac-type phase 5, as one 
can see from (3.54). Hence a proper redefinition of the phases for muon and tau fields 
is needed, in order to make V^ 3 and V t3 real in the course from (3.53) to (3.55). The 
instructive relations between the effective mixing angles in matter 
fundamental mixing angles in vacuum (9 X , 9 y , 9 Z ) are found to be 13 



^Xl 9yi 9 Z 



and the 



tan0 x N 2 + A 




(ml - mf^ 


s\ + (ml - ml) 


2 2" 


D 1 


tan 0s Ny + A 


(ml - m 2 2/ 


si + (ml - m 2 3 ) 


2 2" 

S x C z 


D 2 



tan 9 y 
tan 9 y 

sin 9 Z 
sin#. 



= 1 + 



AlS.m\ x s x c x s z cos 5/ (s y c y 



N 3 -A 

N 3 A 
D~3 + D~ 3 



m\ 



(ml - ml) 



1) 4 + (m 2 3 - ml) c 
+ (m 



m C 



c. 



(3.56) 



12 Note that there is a typing error in Eq. (18) ofRef. [101], associated with the sign of Am| 1 SiciS2e~ 1 ' 5 
in the expression of T1-3. 

13 Here we have only presented the next-to-leading order expression for tan^/tan^, because the 
exact result is rather complicated and less instructive. The former works to a high degree of accuracy, 
and it is identical to the exact result provided that 9 y = 45° and <5 = ±90° hold (i.e., one obtains 



On 



45° in this special but interesting case). 
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Once the relation between 9 y and 9 y is established, one can obtain the relation between 
the effective CP-violating phase 8 in matter and the genuine CP-violating phase 5 in 
vacuum by use of the Toshev identity [102], 

sin (5 = sin 29 y = Sy _ Cy (3 57) 

sin 5 sin 29 y s y c y 

The detailed proof of this elegant identity has been done in Ref. [103]. We see that 5 ~ 5 
holds as a consequence of 9 y ~ 6 y in the leading-order approximation. Hence both 9 y and 
S are insensitive to the matter effects. Note that the results obtained above are only valid 
for neutrinos propagating in matter. As for antineutrinos, the corresponding expressions 
can straightforwardly be obtained from (3.53)-(3.57) through the replacements 5 =>- —5 
and A =>- —A. Such formulas should be very useful for the purpose of recasting the 
fundamental parameters of lepton flavor mixing from the matter-corrected ones [101], 
which can be extracted from a variety of long- and medium-baseline neutrino oscillation 
experiments in the near future. 

To illustrate the dependence of 9 X , 9 y , 9 Z , and 5 on the matter effects, we typically 
take Am^ ps Am s 2 UI1 «5x 10 -5 eV 2 and Am^ w Am 2 atm w 3 x 10 -3 eV 2 . In addition, 
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we adopt 6 X ~ 35°, 6 y rs 40°, # 2 ~ 5°, and 5 ~ ±90°. Such a choice for the values 
of the input parameters is indeed favored by current experimental data on solar and 
atmospheric neutrino oscillations. With the help of (3.56) and (3.57) as well as the 
relevant formulas of m 2 , iVj and A given in section 3.2.1, we are then able to calculate 
the ratios 6 X /6 X , y /0 y , Z /9 Z) and 8/5 as functions of the matter parameter A in the 
typical region 10~ 7 eV 2 < A < 10~ 2 eV 2 , which corresponds to the neutrino beam energy 
E in the region 5 MeV < E < 50 GeV for N e rj 1.5 g/cm 3 [104]. The numerical results 
are shown in Fig. 3.3. We observe that the matter effects on 6 y and S are negligibly 
small for A < 5 x 10~ 2 eV 2 , just as indicated by our analytical results. In contrast, 9 Z is 
significantly modified when A > 10~ 3 eV 2 ; and 9 X is sensitive to the matter effect even 
for quite small values of A. 



3.3 Comparison with Quark Flavor Mixing 



We have seen that the 3x3 lepton flavor mixing matrix V is of a bi-large mixing pattern, 
and most of its matrix elements do not have a strong hierarchy in magnitude. In contrast, 
the 3x3 quark flavor mixing matrix Vckm [10] is of a tri-small mixing pattern and has 
a strongly hierarchical structure. Hence Vckm can be expanded in powers of a small 
parameter A ~ 0.22, the so-called Wolfenstein parameter [105]: 



Vckm 



1-±A 2 
-A 



VAA 3 (1 - p-irj) 



A 

1-|A 2 
-AX 2 




(3.58) 



where the terms of or below C(A 4 ) have been neglected [106]. Current experimental data 
yield A 0.83, p ~ 0.17 and i] ss 0.36 [6]. The spirit of lepton-quark similarity motivates 
us to consider whether the lepton flavor mixing matrix can also be expanded in powers 
of a Wolfenstein-like parameter, whose magnitude may be somehow larger than A. We 
find that it is indeed possible to expand V in powers of A = |V M3 | ~ l/y/2 [91], m view 
of current neutrino oscillation data. 

The parameter A measures the strength of flavor mixing in atmospheric neutrino 
oscillations. The magnitude of flavor mixing in solar neutrino oscillations is then charac- 
terized by AA 2 , where A is a positive coefficient of 0(1). Because |V e3 | < 0.16 is required 
by the CHOOZ data [7], we may take |V~ e3 | ~ 0(A 8 ) ~ 0.06 as a typical possibility for 
A ~ l/y/2 . Smaller values of |V~ e3 | are certainly allowed. In some interesting models of 
lepton flavor mixing [12], \V e3 \ ~ ^Jm e /m^ ~ 0.07 is naturally predicted. Hence C(A 8 ) 
could be the plausible order of \V e ^\. Let us fix the matrix elements V e 2, V e 3 and by 
use of four independent parameters: V^ 3 = A, V e2 = AA 2 and V e3 = BA 8 e~ tS , where the 
positive coefficient B is of 0(1) or smaller and 5 denotes the Dirac phase of leptonic CP 
violation. Then one may make use of the unitarity of V to work out the exact analytical 
expressions for the other six matrix elements. We find that it is more instructive to 
approximate V as 



V = 



A 3 



VI - A 2 A 4 
-AAWl - A 2 
A - BA 5 y/(l - A 2 ) (1 - A 2 A 4 ) e* 



AA 2 



BA H e 



8„-iS 



— A 2 ) (1 - A 2 A 4 ) 



-AVI - ^ 2 A 4 VI - A 2 J 



(3.59) 
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In this approximation, the unitary normalization relations of V keep valid to 0(A ir ) ~ 
2%. Therefore (3.59) is sufficiently accurate to describe lepton flavor mixing, not only 
in solar and atmospheric neutrino oscillations, but also in some of the proposed long- 
baseline neutrino oscillation experiments. As the unitary orthogonality relations of V 
are valid to 0(A 8 ) ~ 6%, the leptonic unitarity triangles and CP violation can also be 
described by (3.59) to an acceptable degree of accuracy. The off-diagonal asymmetries 
of the lepton flavor mixing matrix V [107] read as 





= \Ve2\ 2 - 




Ks\ 2 - 


IK 2 | 2 = 


IK1I 2 - 


\v e3 \ 




= A 2 A 6 , 












A R 


= \Ve2\ 2 - 


K 3 \ 2 = 


Ki\ 2 - 


|K 2 | 2 = 


IKsl 2 - 


\Vel\ 



= A 2 (A 2 A 2 -l) . (3.60) 

We see that Al > holds definitely. In comparison, the sign of Ar cannot be fixed 
from the present experimental data. It is actually possible to obtain Ar = 0, when 
A 2 A 2 = 1 is satisfied. In this interesting case, the lepton flavor mixing matrix V is 
exactly symmetric about its V^-V^-^-i axis [107]. Note that the CKM matrix of quark 
flavor mixing is approximately symmetric about its V u d-V cs -Vtb axis [108]. This is another 
difference between the flavor mixing matrices of leptons and quarks. 

Note that the Majorana CP- violating phases have been omitted in (3.59). To incor- 
porate V with two Majorana phases p and a, we simply multiply V on its right-hand 
side with a pure phase matrix; i.e., V =>- VQ with Q = Diag{e* p , e la , e l5 }. The chosen 
phase convention of Q is to make the Dirac CP-violating phase 5 not to manifest itself 
in the effective mass term of the neutrinoless double beta decay: 

(m) ee = \m 1 (l-A 2 A 4 )e 2ip + m 2 A 2 A 4 e 2ia + m 3 B 2 A 16 \ , (3.61) 

where m; (for % — 1,2, 3) are physical neutrino masses. This result can somehow get 
simplified, if a specific pattern of the neutrino mass spectrum is assumed. The present 
experimental upper bound is (m) ee < 0.35 eV (at the 90% confidence level [52]), from 
which no constraint on p and a can be got. 

Now let us establish the direct relations between (A, A, B) and (0 atm , 6> sun , # chz ). 
With the help of (3.14) and (3.59), we find 

A = sin atm , 
A = 



cos 2 9 chz - \J cos 4 9 chz - sin 2 26» sun 
v^sin 2 ^ 



„ sin 9 C \, 7 

B = — jj-p- . 3.62 

Sin 8 #atm 

Once the mixing angles a tm, 9 sun and 9 C ^ Z are precisely measured, we may use (3.62) to 
determine the magnitudes of A, A and B. For the purpose of illustration, we typically 
take 0.25 < sin 2 sun < 0.40 [35], sin 2 20 atm > 0.92 [2] and sin 2 2# chz < 0.1 [7] to calculate 
the allowed regions of A, A and B. Then we obtain 0.6 < A < 0.8, 0.8 < A < 1.94 
and < B < 9.5. Although the allowed region of B is rather large, we expect that 
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B ~ 0{1) or a bit smaller is most likely and most suitable for our parametrization. It is 
also worthwhile to connect (A, A, B) to (6 X , 6 y , 8 Z ), which are three mixing angles in the 
standard parametrization of V in (3.23). We find sin9 x m AA 2 , sin6 y ~ A and sin #2 = 
BA 8 . In addition, the Dirac phase of CP violation in the standard parametrization is 
exactly equal to 8 defined in the present Wolfenstein-like parametrization. 

An interesting point is that the effective lepton flavor mixing matrix in matter can 
similarly be parametrized in terms of four matter-corrected parameters A, A, B and 8: 



( 



V 



V 7 ! -i 2 A 4 



AA 2 



BA s e~ fs \ 



-iAVl - A 2 



(1 - A 2 )(l -i 2 A 4 ) 



A 



A 3 



I - BA 5 <J(1- A 2 )(l-i 2 A 4 ) e 1 ' 



-AVi - -4 2 A 4 Vi- A 2 



(3.63) 



where the Majorana phase matrix Q has been omitted for simplicity. There also exist 
simple relations between the effective mixing angles of V (i.e., 9 X , 9 y and 6 Z ) and the 
corresponding new parameters (A, A and B). It has been shown in section 3.2.3 that 
sin# y ps smdy and sin 8 ~ sin<5 hold to leading order for a variety of terrestrial long- 
baseline neutrino oscillation experiments. Therefore, we have A A and 8 ~ 8. This 
result implies that A and 8 are essentially stable against terrestrial matter effects. Hence 
the expansion of V in powers of A ps A makes sense. Only A and B in V are sensitive to 
the matter-induced corrections. Because of A oc sin^ and B oc sin# 2 , two remarkable 
conclusions can be drawn from Ref. [103] for our new parameters: (a) A/ A is suppressed 
up to the order Am 2 un / Am 2 tm ; and (b) B/B may have the resonant behavior similar to 
the two-neutrino MSW resonance 



Finally we give some speculation on the physical meaning of A. It is well known that 
the Wolfenstein parameter A ~ 0.22 can be related to the ratios of quark masses in the 
Fritzsch ansatz of quark mass matrices [65] or its modified versions 



A 



i(j>X m d 



(3.64) 



where 0a denotes the phase difference between the (1,2) elements of up- and down- 
type quark mass matrices. This relation indicates that the smallness of A is a natural 
consequence of the strong quark mass hierarchy. Could the largeness of A be attributed to 
a relatively weak hierarchy of three neutrino masses? The answer is indeed affirmative for 
the Fritzsch texture of lepton mass matrices in (2.39), which is compatible with current 
experimental data on neutrino oscillations if the masses of three neutrinos have a normal 
but weak hierarchy (typically, mi : m 2 : m 3 ps 1 : 3 : 10) [66]. In this phenomenological 
model, we approximately obtain 



A 




V rn T 



(3.65) 



where 0a denotes the phase difference between the (2,3) elements of charged lepton and 
neutrino mass matrices. We find that 0a ~ ±180° is practically favored [66], in order 
to obtain a sufficiently large A. To illustrate, we typically take m 2 /m 3 ~ 0.3 as well as 
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m^/rrir ~ 0.06 [6]. Then we arrive at A ~ 0.8, a result consistent with our empirical 
expectation for the order of A in (3.59) 14 . 



4 CP and T Violation 

4.1 Rephasing Invariants in Matter 

No matter whether massive neutrinos are Dirac or Majorana particles, the effects of CP 
or T violation in normal neutrino-neutrino and antineutrino-antineutrino oscillations are 
measured by a universal parameter J in vacuum or J in matter [80], defined through 

lm.{y ai VpjV^Vp i ) = </ ]T e a/37 5> ijfc , 

7 k 

ImiVaiVftVZjVfc) = JE^E 6 ^* (4-1) 

7 k 

where the subscripts (a, /5, 7) and (i,j,k) run respectively over (e,/i, r) and (1, 2, 3). 
Clearly J and J are rephasing-invariant; i.e., they are independent of any redefinition of 
the phases for charged lepton and neutrino fields. Note that J and J can be expressed 
in terms of the moduli of four independent matrix elements of V and V, respectively, as 
follows [12] 15 : 



J 2 = |Ki| 2 |^/|K/|^| 2 -^(i + |Ki| 2 |^| 2 + \v aj \ 2 \Vp. 



. 2 I 2 I 2 1 1 2 ^ 
\Vai\ — \Vf3j\ — \V a j\ — \Vf3i\ 



.I a -l(l + fc 



j 1 = wj'Wai' n;,,r iuf -^n - i^ri^r + i^ri^r 

-|V^| 2 - fef - |K,f - |^| 2 ) 2 , (4.2) 

in which a ^ (3 running over (e, //, r) and 2 7^ j running over (1,2,3). The implication 
of this result is obvious: the information about leptonic CP violation can in principle be 
extracted from the measured moduli of the flavor mixing matrix elements. 

Now let us establish the relationship between J and J. Note that the imaginary parts 
of the rephasing-invariant quantities Z^Z^Z^ and Z e)M Z^ T Z Te , 

lm(Z e ,Z, T Z Te ) = EEE[m 2 m 2 m 2 Im (V^K^IV^)' , 
i=i j=i k=i 

lm(Z e XrZre) = E E E \m 2 m)m\ Im (V^V^V;^] , (4.3) 
i=i j=i k=i 



14 Assuming a somewhat stronger mass hierarchy for three neutrinos, Kaus and Meshkov [109] 
have proposed a different expansion of the neutrino mixing matrix in terms of A = y/m^Jra^ = 
(Aro^/Am^) 1 / 4 ~ 0.37. This parameter is associated with V e 2 instead of V^, therefore it is sensitive 
to the matter effect. In contrast, our paramctrization does not rely on the assumption of neutrino mass 
hierarchy, and its expansion parameter is insensitive to the matter-induced corrections. 

15 Note that there was a typing error in the original formula of J 2 given in Ref. [110]. 
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which do not vanish unless leptonic CP and T are good symmetries, amount to each 
other as a trivial consequence of (3.50). The right-hand side of (4.3) can be expanded in 
terms of J and J. In doing so, one needs to use (4.1) as well as the unitarity conditions 
of V and V frequently. After some lengthy but straightforward algebraic calculations, we 
arrive at an elegant relation between the universal CP-violating parameters J in vacuum 
and J in matter [98] : 

j uffi- ™?) = j n H - m ?) • ( 4 - 4 ) 

i<j i<j 

This interesting result was first obtained by Naumov [111]. Of course J = J holds if 
A = 0, and J = holds if J = 0. It is worth mentioning that (4.4) can also be derived 
from the equality between Det(X) and Det(X) given in (3.46) [112]. Indeed it is easy to 
show that 



2 2 

rri: — m, 
1 j 



Det(X) = 2J J] (m 2 a - m 2 ) J[ 

a</3 i<j 

Det(X) = 2 J (ml - m 2 ) J] (m? - m?) , (4.5) 

a<(3 i<j 

where the Greek indices run over (e,/i, r); and the Latin indices run over (1,2,3). It 
should be noted that the determinant Det(X) or Det(X) contains the same information 
about leptonic CP violation as the universal parameter J or J, although their expressions 
are somehow different [113]. The latter is apparently simpler and more instructive for 
the description of leptonic CP violation in neutrino oscillations. Taking the standard 
parametrization of V in (3.23) or V in (3.55), we have 

J = sin 9 X cos 9 X sin 9 y cos 6 y sin 9 Z cos 2 9 Z sin 5 , 

J = sm9 x cos9 x sm9 y cos9ySm9 z cos 2 9 z sm5. (4.6) 

Thus the CP- violating phases S and S can be related to each other via (4.4) and (4.6). 

Once again, the formulas obtained above are valid only for neutrinos propagating in 
vacuum and in matter. They will become valid for antineutrinos, if the straightforward 
replacements V =^ V* and A =^ —A are made. 

We remark that the effective CP-violating parameter J depends not only upon the 
fundamental CP- violating parameter of the lepton flavor mixing matrix (J), but also 
upon the mass-squared differences of neutrinos (Am^ and Am^) and the matter- induced 
effect (A). In particular, both the sign and the magnitude of J are dependent on the 
signs of Am^ and Am^. With the help of (3.37) and (3.38), we find that the effective 
mass-squared differences Am^ and Am^ can keep unchanged in proper arrangements 
of the signs of Am^, Am^ and A. A careful analysis of (4.4) leads to the following 
exact relations [93]: 

J(+Am 2 1 ,+Am 2 1 ,+A) = - J(-Am 2 21 , -Am 2 31 , -A) , 

J(+Am 2 21 ,-Am 2 31 ,+A) = - J(-Am 2 21 , +Am 2 31 , -A) , 

J(+Am 2 21 ,+Am 2 31 ,-A) = -J(-Am 2 21 ,-Am 2 31 ,+A), 

J(+Am 2 21 ,-Am 2 31 ,-A) = - J(-Am 2 21 , +Am 2 31 , +A) . (4.7) 
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1(T 7 1CT 6 1(T 5 1CT 4 1CT 3 1(T 2 
A (eV 2 ) 

Figure 4.1: The ratio J/ J changing with the matter parameter A for neutrinos (u) and 
antineutrinos (u), in which Am^ = 5 x eV 2 , Am^ = 3 x 10~ 3 eV 2 , 9 X ps 35°, 
6^ 40°, # 2 ~ 5° and 5 ps ±90° have typically been input. 

The validity of these relations are independent of both the neutrino beam energy and 
the baseline length. Note that current experimental data on solar neutrino oscillations 
favor Am^ > 0, but the sign of Am^ remains unknown. 

It is also worth remarking that the CP- and T-violating effects are measured by 
the same parameter (J in vacuum or J in matter), as a straightforward consequence of 
CPT invariance. As for neutrino oscillations, a signal of CP violation comes from the 
probability asymmetry between v a — > vp and V a — > Vp transitions; while a signal of T 
violation is attributed to the probability asymmetry between v a — > vp and vp — > v a 
transitions or between V a — > Vp and Vp — > V a transitions (a ^ (3). The former is 
roughly associated with the difference between J (J, A) and J(—J,—A), in which the 
A-induced effects essentially add each other; but the latter is roughly associated with 
the difference between J(J,A) and J(—J,A) or between J(J,—A) and J(—J, —A), in 
which the A-induced effects essentially cancel each other. For this reason, the CP- 
violating asymmetries in long-baseline neutrino oscillations are rather sensitive to the 
terrestrial matter effect; while the T-violating asymmetries are almost independent of 
the terrestrial matter effect (see Refs. [93, 97] for detailed discussions). Measuring T 
violation is practically more difficult than measuring CP violation, however [114]. 

To illustrate the dependence of J on matter effects, we typically take Am^ = 5 x 
1(T 5 eV 2 and Am 2 31 = 3 x KT 3 eV 2 as well as 6 X ps 35°, B y ^ 40°, 6 Z ps 5° and 
5 ~ ±90° in the standard parametrization of V. The result is shown in Fig. 4.1, where 
1CT 7 eV 2 < A < 1CT 2 eV 2 has been chosen for the matter parameter A. One can see 
that the magnitude of J decreases, when the matter effect becomes significant (e.g., 
A > 1CT 4 eV 2 ). Nevertheless, this feature does not necessarily imply that the CP- 
violating asymmetries in realistic long-baseline neutrino oscillations would be smaller 
than their values in vacuum. Very large terrestrial matter effects can significantly modify 
the frequencies of neutrino oscillations and thus enhance (or suppress) the genuine signals 
of CP violation. 
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Figure 4.2: Unitarity triangles of the lepton flavor mixing matrix in the complex plane. 
Each triangle is named by the index that does not manifest in its three sides [12]. 

4.2 Leptonic Unitarity Triangles 

A geometric description of the lepton flavor mixing phenomenon is physically intuitive 
and instructive. The unitarity of the 3x3 flavor mixing matrix V in vacuum can be 
expressed by two sets of orthogonality relations and two sets of normalization conditions 
for its nine matrix elements: 

i 

E (Ki V *j) = *y, (4-8) 
a 

where Greek and Latin indices run over (e,/i, r) and (1,2,3), respectively. In the com- 
plex plane the six orthogonality relations in (4.8) define six triangles (A e , A M , A T ) and 
(Ai, A 2 , A 3 ) shown in Fig. 4.2, the so-called unitarity triangles. In general, these six 
triangles have eighteen different sides and nine different inner (or outer) angles. The 
unitarity requires that all six triangles have the same area amounting to J/2, where J 
is just the rephasing-invariant measure of CP violation defined in (4.1). If CP were an 
exact symmetry, J = would hold and those unitarity triangles would collapse into 
lines in the complex plane. Note that the shape and area of each unitarity triangle are 
irrelevant to the nature of neutrinos; i.e., they are the same for Dirac and Majorana 
neutrinos. 

In additon to J, there exist other two characteristic quantities of V, resulting from 
its normalization conditions in (4.8). They are the off-diagonal asymmetries Al and Ar. 
defined in (3.60). Clearly Al = (or Ar = 0) would imply that the flavor mixing matrix 
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V were symmetric about its V e i-V^2-V T z (or V^-V^-Ki) axis. Geometrically this would 
correspond to the congruence between two unitarity triangles; i.e., 

A L = =}► A e = Ai , A M = A 2 , A r = A 3 ; 

A* = =S> A e = A 3 , A M = A 2 , A r = A x . (4.9) 

Our current knowledge on lepton flavor mixing is rather poor, hence it remains difficult to 
quantify the magnitude of leptonic CP violation. Nevertheless, we are sure that Ah > 
definitely holds. The possibility for Ar = cannot be ruled out at present [107]. Useful 
information about the sides and angles of the leptonic unitarity triangles will be obtained 
from the medium- and long-baseline neutrino experiments in the near future. 

For instance, one might be able to determine three inner angles of the unitarity 
triangle A 3 , which are defined as 



a x = arg 
Pi = arg 



e2 



V T \V 7 



r2 



* s arg (~W£) ■ (4 - io) 

Then it should be possible to test the self-consistency of our lepton flavor mixing and 
CP-violating picture (e.g., a t +A+7/ 71 would imply the presence of new physics which 
violates the unitarity of the 3x3 lepton flavor mixing matrix). The measurement of a l , 
Pi and r y l can in principle be realized in a long-baseline neutrino oscillation experiment. 
Taking Am^ « Am^ and Am^ « Am^ « ±Amj tm , we simplify (3.26) to [12] 

P{ve^v») = 4|U e3 | 2 |V; 3 | 2 sin 2 F atm - 4Re(v; 1 \/ M2 V;2V; i )sin 2 F sun 



— 8JsinF sllT , sin 2 F. 



atm j 



PW^Vr) = 4|V; 3 | 2 |K 3 | 2 sin 2 F atm - 4Re(\/ Ml K 2 V;* 2 V;*i)sin 2 F sl 



- 8 J sin F sun sin 2 F atm , 
P^-^) = 4|K3| 2 |K 3 | 2 sin 2 P atm - 4Re(KiK 2 V;*2K*i)sm 2 Psun 

- 8JsinP sun sin 2 P atm , (4.11) 

where P atm = P 32 ~ P 3 i and P sun = P 2 i defined below (3.5) measure the oscillation 
frequencies of atmospheric and solar neutrinos, respectively. Provided the baseline 
satisfies the condition L ~ E/Am s 2 mi (i.e., P sun ~ 1), the CP-conserving quantities 
Re(V r Q ,iVa 2 V^ 2 V^ 1 ) and the CP-violating parameter J could both be determined from 
(4.11) for changing values of the neutrino beam energy E. Then one may get 



tana z = 



J 



tan Pi 



Re(V el V^V: 2 V^) ' 
J 



tan 7, = ; . (4.12) 

11 Re(V Tl V e2 V T * 2 V e \) V ' 
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In practice, the earth-induced matter effects are likely to fake the genuine CP-violating 
signals and have to be taken into account. If the neutrino beam energy is low enough 
(e.g., E ~ 0.1 to 1 GeV) and the baseline length is not too long (e.g., L ~ 100 km), the 
terrestrial matter effects might not be very significant for our purpose [115]. 

Another interesting unitarity triangle is A r , whose three sides might more easily 
be measured. To establish A T needs very precise data, which should be able to show 
IKi^iI + IK*3^I > IK*2^ 2 | or \V: 2 V^\ + |K*3^3| > IKi^iI [US]. Such an accuracy 
requirement is practically a big challenge, because one side of A r is much shorter than 
its other two sides (i.e., ~ IK^^I ^ iKIt^l) as (3-19) has shown). Note that 

three sides of the effective unitarity triangle A r in matter, defined as 



A r : Ki^i + V: 2 V, 2 + K* 3 ^3 = , (4.13) 
are possible to be comparable in magnitude. With the help of (3.52), we arrive at 

V ^ ~ Am^Aml/ 61 ^ 1 + Am 2 2X Am 2 J^ 2 ' 



~ ~ _ Am 2 2 Am 2 21 Am 2 22 Am 2 31 

V ^ 2 ~ Am 2 2 Am 2 J^ 2 + Am 2 2 Am 2 21 V « V >* ' 

~ ~ Am 2 3 Am 2 23 Am 2 33 Am 2 21 

Ka ^ 3 ~ AmlArh 2 /^ 3 + A^Amf/^ 1 ' (4 " 14) 

where Am 2 - = m 2 — m 2 . The explicit expressions of Am 2 - and Am 2 -, which depend on 
both Am 2 - and A, can be obtained from (3.37). The instructive result in (4.14) clearly 
shows how three sides of A T get modified in matter. Of course, one may use (4.14) to 
evaluate three sides of A T . It is also possible to get a numerical feeling of A r and A T 
directly from Fig. 3.2. For example, we find that three sides of A r become comparable 
in magnitude (of order |V e * 3 V^3|), when the neutrino beam energy E is about 1 GeV or 
somehow larger. In a similar way, we can define and discuss other effective unitarity 
triangles in matter. 

If \V e3 \, \Vfja\ and \V t3 \ are well determined in the first-round long-baseline neutrino 
experiments, it should be possible to check one of the three normalization conditions: 
| V^3 1 2 + |V^ 3 | 2 + IK3I 2 = 1- If \V e i\, |V e2 | and \V e s\ are measured to a good degree 
of accuracy, one can test another normalization condition \V e i\ 2 + | V^2 1 2 + |Ve3| 2 = 1- 
At this stage the moduli of V^i, V^, V T \ and V T 2 remain unknown. To determine the 
universal CP-violating parameter J and to make a full test of the unitarity of V, much 
more delicate long-baseline neutrino experiments are needed. 



4.3 Baryon Asymmetry via Leptogenesis 

In the universe, the density of baryons compared to that of photons is extremely small: 
7] = n B /n 7 = (2.6 — 6.3) x 10~ 10 , extracted from the Big-Bang nucleosynthesis [6]. This 
tiny quantity measures the observed matter-antimatter or baryon-antibaryon asymmetry 
of the universe, 

^ n B ~n^ ^ rj_ = _ ^ 1Q _ n 
s 7.04 
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where s denotes the entropy density. To produce a net baryon asymmetry in the standard 
Big-Bang model, three Sakharov necessary conditions have to be satisfied [117]: (a) 
baryon number nonconservation, (b) C and CP violation, and (c) a departure from 
thermal equilibrium. Among a number of interesting and viable baryogenesis mechanisms 
proposed in the literature [118], Fukugita and Yanagida's leptogenesis mechanism [119] 
has recently attracted a lot of attention - due partly to the fact that neutrino physics is 
entering a flourishing era. 

As mentioned in section 2.1.2, a simple extension of the standard model to gen- 
erate neutrino masses is to include one right-handed neutrino in each of three lepton 
families, while the Lagrangian of electroweak interactions keeps invariant under the 
SU(2) L x U(l)y gauge transformation. In this case, the Yukawa interactions of leptons 
are described by 16 



-C Y = lh<j>Yie R + l L (j)Y u v R + -v R M R v R + h.c. 



(4.16) 



where /l denotes the left-handed lepton doublet, eR and z/r stand respectively for the 
right-handed charged lepton and Majorana neutrino singlets, and is the Higgs-boson 
weak isodoublet. The lepton number violation induced by the third term of Cy allows 
decays of the heavy (right-handed) Majorana neutrinos iVj (for i = 1,2, 3) to happen; 
i.e., Ni — > I + 0* vs Ni — > l c + 0, as illustrated in Fig. 4.3. Because each decay mode 
occurs both at the tree level and at the one-loop level (via the self-energy and vertex 
corrections), the interference between the tree-level and one-loop decay amplitudes can 
lead to a CP-violating asymmetry £j between the two CP-conjugated processes [121]: 



Si = 



i + 0t) _ r(Ni -> i c + 0) 



T(Ni -> I + 0t) + r(Ni - l c + 0) 



8vr MyJ) . 



V 



'Mf 



fs 




(4.17) 



where the loop functions / v (vertex) and / s (self-energy) are given by 

1 +x] 



fv(x) = 



fs(x) = < 



If x ^> 1 holds, we arrive at / s 



1 - (1 + x)ln 
1 + x 



\fx In 

x 



1 - x 

2y^ 

1 -x 
2/ v 



X 

(SM), 



x 

(SUSY) 



(SM), 



(SUSY) . 



2/ 



v 



(4.18) 
-2/y/E (SUSY). 



-1/v^ (SM) or f s 

The central idea of baryogenesis via leptogenesis is rather simple [119]. In the early 
universe, the heavy (right-handed) Majorana neutrinos iVj are present in the primordial 

16 One may also consider a similar extension of the minimal supersymmetric standard model to gen- 
erate neutrino masses and to resolve the hierarchy problem induced by heavy right-handed Majorana 
neutrinos [120]. The parameter counting for lepton flavor mixing and leptogenesis is identical in models 
with and without supersymmetry [24]. 
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(a) 





Figure 4.3: Feynman diagrams for the decay of a heavy (right-handed) Majorana neutrino 
Ni\ (a) the tree-level contribution; (b) the one- loop self-energy correction; and (c) the 
one-loop vertex correction. Note that JVj's antiparticle is itself up to a phase factor, and 
"x" indicates the chirality flip between JVj and its antiparticle in the propagator. For 
simplicity, we have only taken account of the diagrams which may result in the CP-odd 
contributions to the decay rates of JVj. 



thermal bath. They may decay into leptons and scalars as soon as the universe becomes 
cold enough. As shown in Fig. 4.3, the decays of JVj violate both the lepton number 
conservation and the CP symmetry. If such decays occur out of thermal equilibrium, a 
net lepton-antilepton asymmetry can be generated. Later on, this net lepton-antilepton 
asymmetry is converted in part to a net baryon-antibaryon asymmetry through the 
(£>+L)-violating sphaleron processes [122]. For simplicity, we consider the case that three 
heavy Majorana neutrinos JVj have a hierarchical mass spectrum (Mi <C M2 <C M3). 
Then the interactions of JVi are in thermal equilibrium when JV2 and JV3 decay, and the 
CP-violating asymmetries produced in the decays of JV 2 and JV 3 (i.e., £2 and e 3 ) can 
be erased before JVi decays. The CP-violating asymmetry S\ produced by the out-of- 
equilibrium decay of JVi survives, and it results in a net lepton-antilepton asymmetry 

Y L s = l £l , (4.19) 

where g # = 106.75 (SM) or 228.75 (SUSY) is an effective number characterizing the 
relativistic degrees of freedom which contribute to the entropy s, and d accounts for 
the dilution effects induced by the lepton-number- violating wash-out processes [121]. 
The lepton asymmetry Y L is eventually converted into a net baryon asymmetry Y B via 
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nonperturbative sphaleron processes [122]. The explicit relation between Yl (initial) and 
Yb (equilibrium) is given as [123] 



where c = (8N/ + 4N^)/(22N/ + 13N^) with Nj being the number of fermion families 
and being the number of Higgs doublets. Taking N/ = 3 and = 1 for example, 
we obtain c ~ 1/3. Note that the dilution factor d can be computed by integrating 
the full set of Boltzmann equations [124]. In the literature [125], some useful analytical 
approximations for d have frequently been made. 

In the simple mechanism of thermal leptogenesis introduced above, the baryon asym- 
metry of the universe is attributed to the out-of-equilibrium decay of the lightest heavy 
(right-handed) Majorana neutrino N\. Another interesting scenario of thermal leptogen- 
esis is to consider the lepton- number- violating decays of two heavy Majorana neutrinos, 
whose masses are approximately degenerate [126]. Because the self-energy contribution 
to £j can significantly be enhanced in the case of Mj m Mj, it is possible to generate 
the observed baryon-antibaryon asymmetry Yb via the out-of-equilibrium decays of two 
relatively light and moderately degenerate iVj. Such a scenario could allow the smallest 
mass of Ni to be close to or below the maximum reheating temperature of the universe 
after inflation in the generic supergravity models. Of course, one may abandon the idea 
of thermal leptogenesis and interpret the baryon asymmetry through the decays of heavy 
(right-handed) Majorana neutrinos produced non-thermally by the inflaton decay [127]. 
The non-thermal leptogenesis mechanism seems to be more speculative and less elegant 
than the thermal leptogenesis mechanism, however. Thus we continue to focus on the 
simple scenario of thermal leptogenesis in the following. 

After spontaneous symmetry breaking, C Y in (4.16) becomes C mass in (2.6), where 
Mi = Yi((f>) for the charged lepton mass matrix and Mq = Y v ((f>) for the Dirac neutrino 
mass matrix. The scale of M/ and Md is characterized by the gauge symmetry break- 
ing scale v = (0) ~ 174 GeV; but the scale of Mr may be much higher than v, because 
right-handed neutrinos are SU(2) L singlets and their mass term is not subject to the elec- 
troweak symmetry breaking. As already shown in section 2.1.2, the light (left-handed) 
neutrino mass matrix M v can be given in terms of M^ and Mr via the seesaw relation 
(2.9). Note that lepton flavor mixing at low energy scales stems from a nontrivial mis- 
match between the diagonalizations of M v and M h while the baryon asymmetry at high 
energy scales depends on complex Y v and Mr in the thermal leptogenesis mechanism. 
To see the latter more clearly, let us diagonalize the symmetric mass matrix Mr by a 
unitary transformation: 



where Mi are the physical masses of JVj. Provided three heavy Majorana neutrinos have 
a strong mass hierarchy (i.e., M x <C M 2 <C M 3 ), the CP-violating asymmetry between 
iVi -> / + 1 " and Nx -> l c + <\> decays in (4.17) can be simplified to [128] 



Y B = -cY L = -c—Ei , 

9* 



(4.20) 




(4.21) 



16ttv 2 [[/^M^M D C/ R ] 




ii 




(4.22) 
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In obtaining (4.22), we have taken into account the point that the Dirac neutrino Yukawa 
coupling matrix Y u takes the form MdII^/v in the physical basis where Mr is diagonal. 
One can see that S\ depends on the complex phases of Md and C/r. The only possible 
relationship between E\ and the lepton flavor mixing matrix V at low energy scales is due 
to the seesaw relation (2.9), which links M v to Md and Mr. Therefore we conclude that 
there is no direct connection between CP violation in heavy Majorana neutrino decays 
[ei) and that in neutrino oscillations (J). Such a general conclusion was first drawn by 
Buchmuller and Plumacher in Ref. [129]. Recently some other authors have carried out 
more delicate analyses and reached the same conclusion. 

Depending on the specific flavor basis that we choose in model building, and V 
can either be completely disconnected or somehow connected. To illustrate, we consider 
two extreme cases [130]: 

• In the flavor basis where M v is diagonal (i.e., lepton flavor mixing and CP violation 
at low energy scales arise solely from the charged lepton sector [61]), we find that 
E\ has nothing to do with V. In this special case, less fine-tuning is expected in 
building a phenomenological model which can simultaneously interpret the baryon 
asymmetry of the universe and lepton flavor mixing at low energy scales. 

• In the flavor basis where both M/ and Mr are diagonal (i.e., lepton flavor mixing 
and CP violation at low energy scales arise solely from the neutrino sector), we 
find that e\ can indirectly be linked to V through the seesaw relation in (2.9). In 
this special case, it is highly nontrivial to build a predictive model or ansatz which 
can simultaneously interpret the observed baryon asymmetry of the universe and 
current neutrino oscillation data. 

It is worth remarking that the correlation between high-energy physics and low-energy 
physics may offer us a valuable opportunity to probe the former from the latter, or vice 
versa [131]. When a specific model is built, however, the textures of Mr and Md have to 
be carefully chosen or fine-tuned to guarantee acceptable agreement between the model 
predictions and the observational or experimental data (see Refs. [128, 132] for example). 
Many different models or ansatze on leptogenesis and neutrino oscillations have recently 
been proposed [133], but some of them are rather preliminary and speculative. 

An important question is whether the leptogenesis mechanism can experimentally be 
proved. The answer is unfortunately negative. But baryogenesis via leptogenesis might 
become a conceivable and acceptable mechanism to interpret the matter-antimatter 
asymmetry of the universe [134], if (1) the electroweak baryogenesis scenario is ruled 
out; (2) the Majorana nature of massive neutrinos is verified (e.g., by measuring the 
neutrinoless double beta decay); and (3) the leptonic CP- violating effect is observed 
(e.g., in the future long-baseline neutrino oscillation experiments). 

5 Concluding Remarks 

We have presented an overview of recent progress in the phenomenological study of 
neutrino masses, lepton flavor mixing and CP violation. Particular attention has been 
paid to the model-independent properties of massive neutrinos in vacuum and in matter. 
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With the help of current experimental and observational data, we have obtained some 
enlightening information about the neutrino mass spectrum. To pin down the absolute 
mass scale of three active neutrinos, one has to make much more efforts to detect the 
tritium beta decay and the neutrinoless double beta decay. Precise cosmological data are 
also expected to play an important role in determining the absolute values of neutrino 
masses. The relative magnitudes of three neutrino masses can be fixed by means of a 
variety of long-baseline neutrino oscillation experiments in the coming years. 

Different from quark flavor mixing, lepton flavor mixing involves two large mixing an- 
gles. Whether the largeness of neutrino mixing angles is associated with a relativly weak 
hierarchy of neutrino masses remains an open question. We hope that both the small- 
est angle of lepton flavor mixing and the Dirac phase of CP violation can be measured 
in the long-baseline neutrino oscillation experiments. It seems hopeless to separately 
determine two Majorana phases of CP violation from the measurements of the neutri- 
noless double beta decay and other lepton-number-violating processes. Nevertheless, it 
is extremely important to realize such measurements in order to identify the Majorana 
nature of massive neutrinos and to shed light on the nontrivial features of lepton number 
violation. 

To conclude, the robust experimental evidence for neutrino masses and lepton flavor 
mixing strongly indicates that the standard electroweak model is actually incomplete. 
This incompleteness motivates us to open a new window to go beyond the standard 
model. Although there have been a number of unresolved questions associated with 
massive neutrinos, we are certainly paving the way for satisfactory answers to them. A 
convincing and predictive theory of massive neutrinos should be achieved in the foresee- 
able future. Such a theory, which must include new physics of both leptons and quarks at 
high energy scales, is expected to offer us a deeper insight into the generation of fermion 
masses, the pattern of flavor mixing and the origin of CP violation. 
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A Mixing between Active and Sterile Neutrinos 

In addition to the robust evidence for atmospheric and solar neutrino oscillations accu- 
mulated from the SK [2], SNO [3], KamLAND [4] and K2K [5] experiments, — > V e and 
v n v e transitions have been observed by the LSND Collaboration [8]. The LSND data 
can also be interpreted in the assumption of neutrino oscillations, whose mass-squared 
difference and mixing factor read as 

AKsnd ~ 1 eV 2 , sin 2 2# LSND ~ 1(T 3 - 1(T 2 . (A.l) 

Because solar, atmospheric and LSND neutrino oscillations involve three distinct mass- 
squared differences (Am 2 un <C Am 2 tm <C Am 2 SND ), a simultaneous interpretation of 
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Figure A.l: Sample mass spectra of two four-neutrino mixing schemes, where each line 
denotes mf (for % — 0, 1, 2, 3). 

them requires the introduction of a light sterile neutrino 17 . There are two typical 
categories of four-neutrino mixing schemes and their sample mass spectra are shown in 
Fig. A.l: (a) the (2+2) scheme with two pairs of nearly degenerate massive neutrinos, 
whose mass-squared gap is characterized by AmL SND ; and (b) the (3+1) scheme with a 
triplet of nearly degenerate massive neutrinos and an isolated massive neutrino, whose 
mass-squared gap is also characterized by AmL SND . Although each four- neutrino mixing 
scheme consists of three independent mass-squared differences and six flavor mixing 
angles, its parameter space can be tightly constrained by current experimental data. The 
global analyses of all available neutrino oscillation data [137] have shown that the (2+2) 
mixing scenario is strongly disfavored, while the (3+1) mixing scenario is marginally 
allowed. However, the recent cosmological upper bound on the sum of all light neutrino 
masses in (2.22) is too low to be consistent with the LSND result given in (A.l), at least 
at the 95% confidence level [48]. It is argued by Hannestad [47] and Giunti [48] that the 
(3+1) mixing scheme may still survive in a very tiny parameter space, if the relevant 
cosmological data are considered at a higher confidence level. The upcoming MiniBooNE 
experiment [39] is crucial to firmly confirm or disprove the LSND measurement. 

Before a definitely negative conclusion can be drawn from MiniBooNE, we do not 
think that the LSND data should be completely discarded. In particular, it is worthwhile 
to study the four-neutrino mixing scenarios in a way without special theoretical biases 
and (or) empirical assumptions. Starting from such a point of view, we present some 
model-independent results on the description of four-neutrino mixing and CP violation 
in this Appendix. 

A.l Standard Parametrization 

Let us begin with a generic SU(2)l x U(1)y model of electroweak interactions, in which 
there exist n charged leptons belonging to isodoublets, n active neutrinos belonging to 
isodoublets, and n' sterile neutrinos belonging to isosinglets. The charged-current weak 
interactions of leptons are then associated with a rectangular flavor mixing matrix of n 

17 A sterile particle means that it has vanishing or extremely feeble couplings to the other particles in 
the standard model. Instead of introducing a light sterile neutrino, a few more far-fetched ideas (such 
as the violation of CPT symmetry in the neutrino sector [135] and the lepton-number-violating muon 
decay [136]) have been proposed in the literature to simultaneously interpret the solar, atmospheric and 
LSND anomalies within the three-neutrino mixing scheme. 
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rows and (n + n') columns [72]. Without loss of generality, one may choose to identify the 
flavor eigenstates of charged leptons with their mass eigenstates. In this specific basis, 
the n x (n + n') lepton mixing matrix links the neutrino flavor eigenstates directly to the 
neutrino mass eigenstates. Although sterile neutrinos do not participate in normal weak 
interactions, they may oscillate among themselves and with active neutrinos. Once the 
latter is concerned we are led to a more general (n + n') x {n + n') lepton flavor mixing 
matrix [12], defined as V in the chosen flavor basis. For the mixing of one sterile neutrino 



(u s ) and three active neutrinos 



ye, Vp, 



v T \ the explicit form of V can be written as 
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where z/j (for i — 0, 1, 2, 3) denote the mass eigenstates of four neutrinos. If neutrinos 
are Dirac particles, V can be parametrized in terms of six mixing angles and three phase 
angles. If neutrinos are Majorana particles, however, three additional phase angles are 
required to get a full parametrization of V. We totally need six complex rotation matri- 
ces, denoted as i?oi, -R02, R03, R12, R13 and R23, which correspond to simple rotations 
with angles % in the (0,1), (0,2), (0,3), (1,2), (1,3) and (2,3) planes. For simplicity, 
we assume that each Rij involves only a single phase angle Sij associated with its sin % 
term. Explicitly, we have 
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where Sij = sin e tSij and = cos % . There exist numerous different ways to arrange 
the products of these rotation matrices [78] , such that the resultant flavor mixing matrix 
V covers the whole 4x4 space. After a proper assignment of the relevant rotation angles, 
one can find that only sixteen parametrizations of V are structurally distinct: four of 
them have the cos 9ij terms in the (i, i) positions of V; and twelve of them have the sin 
terms in the twelve different (i, j) positions of V (for i ^ j). 

To avoid unnecessary complication, we shall not write out a lengthy list of the sixteen 
parametrizations of V. Instead we only take an instructive example for illustration. 
Similar to the representation in (3.23), the 4x4 neutrino mixing matrix V can be 
parametrized as V = R23 <E> -R13 <8> -R03 <8> R12 <8> -R02 <S> R01 [138]; i.e., 



V = 
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Without loss of generality, the six mixing angles ^ can all be arranged to lie in the 
first quadrant. The six CP- violating phases 5^ may take arbitrary values between and 
2tt. Note that V can also be decomposed into a Dirac-like flavor mixing matrix with six 
rotation angles and three CP-violating phases, multiplied by a diagonal phase matrix 
with three Majorana phases. In normal neutrino-neutrino and antineutrino-antineutrino 
oscillations, CP-violating phases of the Majorana nature are unable to be measured. 



A. 2 Invariants of CP Violation 

No matter whether neutrinos are Dirac or Majorana particles, one may define the Jarl- 
skog rephasing invariants of four- neutrino mixing [80], which govern the CP- and T- 
violating effects in normal neutrino-neutrino and antineutrino-antineutrino oscillations. 
To be explicit, we have 

J% = \m(v m V P] V: 3 V^ , (A.5) 

where the Greek subscripts run over (s, e, /i, r) and the Latin superscripts run over 
(0, 1, 2, 3). Of course, J% = = and J% = -J% = -J% = jg a hold by definition. 
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The unitarity of V leads to the following correlation equations of J^: 

T,J% = H J % = Y, J % = T,J% = o- (a.6) 

i j a p 

Hence there are totally nine independent J^, whose magnitudes depend only upon three 
of the six CP-violating phases or their combinations in a specific parametrization of V. 

We use the standard parametrization of V in (A. 4) to calculate J l Jp. After some 
lengthy but straightforward calculations, we obtain the exact analytical expressions of 
nine independent [140]: 

t02 2 2 i 2 2 2 j. 

J TS = C 01 C02C 03 Cl2Cl3C23S02S03S23Sm(p x + CqiCq2C 03 Ci2Ci 3 C 23 SoiS 02 So 3 Si 3 Sin <p y 

+ ( c 23 s i3 _ s 23) c oiCo 2 Co 3 Ci2SoiSo2Si2 sin (f) z - c 01 cl 2 cl 3 cl 2 c 23 s is 02 si 3 s 2 3 sin(^ x - 4> y ) 

-C iCo2Co 3 Ci3C 2 3SoiSo 2 So3Sl2S23 sin(0 x + (j) z ) + Co 1 Co2Co 3 Ci3C 23 So2So3Si2Si3 sin(</> y - <j) z ) 

-coiCo 2 Co 3 c 23 s iSo2Si 2 si3S23 sin(0 x - 4> v + 2(j) z ) , 

03 2 2 2 2 

J TS = -c ol co2C O3 ci2Ci 3 C23So2So3S23sin0 2; - c ic 2C 03 ci2C 13 c 23 soiS 03 Si 3 sin 4> y 

+c ic 2Co 3 ci 3 C23SoiSo3Si2S 2 3 sin((j) x + 4> z ) - C01C02C03C13C23S02S03S12S13 sm((f) y - <j> z ) , 

Jrs = C02Cg 3 Ci2Ci 3 C23S02S03S23 Sin^ + C 2Co 3 Ci3C 23 So2S03Sl2Sl3 Sin((/) y - Z ) ; (A. 7) 

and 

Jse = c OlC02Co 3 Cl2Ci3SoiSo2 s 03Sl3 Sm<j)y - CqiC^ 2 C^ 3 Ci 2 c\ 3 SqiSq2S\2 SHI <j) z 

+C01C02C03C13S02S03S12S13 sm((j>y - <j) Z ) , 
Jll = C lC 2Co 3 Ci2Ci 3 SoiSo3Sl3 Sm <j) y - C 02 Co 3 Ci 3 SoiS02S03Sl2Sl3 sin((/)y - (j) z ) , 
Jse = C O 2Co 3 Ci 3 S O 2SO3Sl2Sl3Sin(0 y - <j) z ) ; (A. 8) 

as well as 

Jll = - (coiCi 2 s? 3 - CoiC? 3 sf 2 - Soi s 03 s i3 + s oi s i2) C02C12C13C23S02S03S23 sin <j> x 
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+ CoiCo2Cl3C23SoiSo 2 S03S? 2 Sl3S23 Sm{(j) x - 2(j) y + 3<f> z ) , 
Jll = C 02 Co 3 Ci2Ci 3 C 2 3SoiSo2S03Si3S23 Sm 4> x + CoiC 02 Co 3 Ci2Ci 3 SoiS03Sl3S23 sin (j) y 

- c OlCo 3 Ci 2 Ci 3 C23SoiS02Sl3S23 S\n{<j) x - <j)y) + CoiCo2Co 3 Ci 3 C23SoiSo3Sl2Si 3 S23 sin((f) x + (f> z ) 
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+ c 02 c 03 c 12 c 13 c 23Sl 2 Sl 3 S23 Sm((j) x - <j> y + <f> z ) , (A.9) 
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where 



4>x = ^03 - <5o2 - #23 , 
4>y = ^03 - ^01 - #13 , 

<Pz = 5 02 - 5 01 - 5 12 . (A. 10) 

With the help of (A. 6), one may easily derive the expressions of all the other rephasing 
invariants of CP and T violation from (A. 7), (A. 8) and (A. 9). The results obtained above 
are expected to be very useful for a systematic study of CP- and T-violating effects in 
the four-neutrino mixing models. The same results are also applicable for the discussion 
of CP and T violation in the four-quark mixing models [77, 139]. 

Note that all CP- and T-violating observables in neutrino oscillations must be related 
linearly to J^p- To see this point more clearly, we consider that a neutrino u a converts 
to another neutrino Vp in vacuum. The probability of this conversion is given by 

P{u a -+ vp) = <W - 4 ]T [Re (VaVpjV^jV^ sin 2 F jt ] - 2 £ ( J% sin 2F jt ) , (A.ll) 

i<j i<j 

where Fji = 1.27 Am^L/E with Amjj = m| — m 2 , L stands for the baseline length 
(in unit of km), and E is the neutrino beam energy (in unit of GeV). CPT invariance 
assures that the transition probabilities P{yp — > u a ) and P(v a — > Vp) are identical, and 
they can directly be read off from (A.ll) through the replacement J^p ==>- —J^p (i.e., 
V =>- V*). Thus the CP- violating asymmetry between P(u a — > Vp) and P(V a — > Vp) is 
equal to the T-violating asymmetry between P(u a — > Vp) and P{yp — > v a ). The latter 
can be explicitly and compactly expressed as follows [140]: 



AP a p = P{vp -> i/J - P{y a -> ^) 

= 16 (j^p sin F 2 i sin F 31 sin F 32 + sin F w sin F 30 sin F 31 

+ J° 2 p sin F 20 sin F 30 sin F 32 ) . (A.12) 
Equivalent ly, one may obtain 

AP a p = 16 ( Jlp sin F 21 sin F 31 sin F 32 - Jfp sin F 10 sin F 20 sin F 21 



J°| sin Fio sin F 30 sin F 31 ) , (A. 13) 



or 



AP a p = 16 ( Jfp sin F 21 sin F 3 i sin F 32 + J°J sin Fi sin F 20 sin F 2i 

- sin F 20 sin F 30 sin F 32 ) . (A. 14) 

In getting Eqs. (A.12)-(A.14), the equality 

sin 2Fij + sin 2F jk + sin 2F ki = — 4 sin sin sin F ki (A. 15) 

and (A. 6) have been used. Only three of the twelve asymmetries AP a p are independent, 
and they probe three of the six CP-violating phases (or their combinations) of V. Since 
only the transitions between active neutrinos can in practice be measured, the useful 
probability asymmetries for the study of leptonic CP and T violation in long-baseline 
neutrino oscillation experiments are AP efl , AP^ T and AP re . 
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A. 3 Leptonic Unitarity Quadrangles 



The unitarity of V implies that there exist twelve orthogonality relations and eight 
normalization conditions among its sixteen matrix elements. The former corresponds 
to twelve quadrangles in the complex plane, the so-called unitarity quadrangles. To be 
explicit, let us write out the twelve orthogonality relations and name their corresponding 
quadrangles [141]: 
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(A.16) 
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If six mixing angles and six CP-violating phases of V are all known, one can plot twelve 
unitarity quadrangles without ambiguities. Note, however, that each quadrangle has 
three distinct topologies in the complex plane. For illustration, we take quadrangle Q se 
for example and show its three topologies in Fig. A. 2, where the sizes and phases of V^l/* 
(for % = 0,1,2,3) have been fixed. One can see that different topologies of quadrangle Q se 
arise from different orderings of its four sides, and their areas are apparently different 
from one another. As a whole, there are totally thirty-six different topologies among 
twelve unitarity quadrangles. 

Now we calculate the areas of all unitarity triangles and relate them to the rephasing 
invariants of CP violation J^. The results for thirty-six different topologies of twelve 
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unitarity quadrangles are summarized as follows 18 : 
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(A.18) 



where the subscripts a/3 = se, s/i, sr, e/i, er or /rr, and ij = 01, 02, 03, 12, 13 or 
23. Note that the correlation of J^p allows us to simplify (A.18) 



Then each S^p or 



(for q = a,b, c) can be expressed as a sum of two independent Jarlskog invariants. Such 
simplified expressions of S^g and Sfj depend on the choice of independent J^ g , therefore 
they may have many different forms. If nine independent Jarlskog invariants are fixed, 
however, some expressions of S^ g and Sfj must consist of three J^p- This point will 
become clear later on. 

As Jap = ~Jap = ~Jp a = Jp a n °lds by definition, one may easily obtain S q a3 = —Sp a 
and Sfj = —S^ (for q = a, b, c). From the sum rule in (A. 6), one can also find 
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where a or (3 runs over (s, e, /i, r), and i or j runs over (0, 1, 2, 3). In addition to (A.19), 
the following relations can be derived from (A.18): 
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18 It should be noted that the areas of unitarity quadrangles under discussion are "algebraic areas" , 
namely, they can be either positive or negative. Of course, it is always possible to take Sg e = (| \ + 
\Jfe\ + \Jfe\ + or Sf e = |Ji° + Jf e + Jf e + J s ° e 3 |/4, such that S« e is dehnitely positive. We 

find, however, that the language of "algebraic areas" is simpler and more convenient in the description 
of unitarty quadrangles. In particular, the algebraic area of unitarity quadrangle Q se in the case of 
topology (b) means an algebraic sum of the areas of its two disassociated triangles, which have opposite 
signs. Hence both Sg e = and S|? e < are in general allowed. As for topologies (a) and (c) of Q se , 
Sg e > and S c se > are simply a matter of sign or phase convention. 
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The correlative relations in (A. 19) and (A. 20) indicate that there are only nine indepen- 
dent S^p and (or) Sfj, corresponding to nine independent J^g- 

Without loss of generality, let us choose the following nine independent S^g [141]: 
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In terms of six flavor mixing angles and three independent phase combinations of V, we 
have expressed the nine independent J^ appearing on the right-hand side of (A.21) 
in section A. 2. Then one may directly obtain the explicit expressions of the nine- 
independent S^g in terms of the same mixing angles and CP-violating phases. 



A. 4 Matter Effects on T Violation 



The effective Hamiltonians responsible for the propagation of active and sterile neutrinos 
in vacuum and in matter can respectively be written as [33] 
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where D v = Diag{mo, mi, m2, m^} with m^ being the fundamental neutrino masses in 
vacuum, D u = Diag{m , mi, m 2 , m 3 } with m, being the effective neutrino masses in 
matter, and E ^> rrii denotes the neutrino beam energy. The deviation of 7i e s from 7i eff 
is given by 
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(A.23) 



where a = a/2 GpN e and a' = \/2 GpN n /2 with iV e and N n being the background 
densities of electrons and neutrons [142], respectively. In the literature one often assumes 
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a constant earth density profile (i.e., N e = constant and N n = constant), which is a good 
approximation for all of the presently-proposed long-baseline neutrino experiments. 

Now let us introduce the commutators of 4 x 4 lepton mass matrices to describe 
the mixing of one sterile and three active neutrinos. Without loss of any generality, 
we continue to work in the afore-chosen flavor basis, where Mj takes the diagonal form 
Di = Diag{m s , m e , m M , m T } with m s = 0. Note that we have assumed the (1, 1) element 
of Di to be zero, because there is no counterpart of the sterile neutrino v s in the charged 
lepton sector. We shall see later on that our physical results are completely independent 
of m s , no matter what value it may take. The commutator of lepton mass matrices in 
vacuum and that in matter can then be defined as 
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Obviously C and C are traceless Hermitian matrices. In terms of neutrino masses and 
flavor mixing matrix elements, we obtain the explicit expressions of C and C as follows: 
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c = % 



where A a p = — m| for a 7^ (3 running over (s, e, /x, r), and 



(A.25) 
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(A.26) 



One can see that A f 
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related to zT Q( g, we need to establish the relation between C and C. Taking account of 
(A. 22) and (A. 23), we immediately obtain 
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This interesting result indicates that the commutator of lepton mass matrices in vacuum 
is invariant under terrestrial matter effects. As a straightforward consequence of C — C, 
we arrive at Z a p = Z a p from (A.25); i.e., 



(A.28) 



i=0 
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or equivalently 



(A.29) 



i=i 



i=l 



where A i0 = to? — m 2 , and Aj = mf — m 2 , for i — 1,2, 3. It becomes obvious that the 
validity of (A. 28) or (A.29) has nothing to do with the assumption of m s = in the 
charged lepton sector. Note that the results obtained above are only valid for neutrinos 
propagating in vacuum and in matter. As for antineutrinos, the corresponding sum rule 
can straightforwardly be written out from (A. 28) or (A.29) through the replacements 
V =^ V*, a =^ —a and a' =^ —a'. 

To describe CP or T violation in neutrino oscillations, we consider the rephasing- 
invariant relationship Z a ^Z Pl Z ia = Z a ^Z^Z ia for a ^ f3 ^ 7 running over (s,e,/i, r). 
The imaginary parts of Z aj gZ J g 7 Z 7Q , and Z a pZ^Z^ a read explicitly as 

i=\ j=i k=i 

lm(Z a pZ^Z, a ) = ]T ]T £ [AioAj-oAw Im (K^ 7 ^4^)], (A.30) 

i=i j=i fc=i 

and their equality allows us to derive an interesting relation between the Jarlskog invari- 
ant in vacuum and its effective counterpart in matter, 



J% = Im (VaiVpjV:^) , 



(A.31) 



where the Greek subscripts run over (s, e, /i, r) and the Latin superscripts run over 
(0,1,2,3). The result is [98] 
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If one "switches off" the mass of the sterile neutrino and its mixing with active neutrinos 
(i.e., a' = 0, A i0 = to 2 , A i0 = mf, = 0, and = 0), then (A.32) is simplified to 
the elegant Naumov form [111], as shown in (4.4). 

The matter-corrected CP-violating parameters J^g can, at least in principle, be deter- 
mined from the measurement of CP- and T- violating effects in a variety of long-baseline 
neutrino oscillation experiments. The conversion probability of a neutrino v a to another 
neutrino v* in matter is given as 



i<j 
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where Fji = \.TJ A^Lj E with Aji = — mf. The transition probability P{y a — > f Q ) 
can directly be read off from (A. 33), if the replacements J^a ==>- —J^b are m &de 19 . To 
obtain the probability P(v a — > 77g), however, both the replacements ==>- —J^ g and 
(a, a') ==>- (—a,— a') need be made for (A. 33). In this case, P(v a — > V g ) is not equal 

to P{y B — > z/q,). The difference between P(z7 a — * Vq) and P(z^3 — > ^«) is a false signal 
of CPT violation, induced actually by the matter effect [73]. Thus the CP- violating 
asymmetry between P{y a — > i/g) and P(z7 a — > V : /3 ) is in general different from the T- 
violating asymmetry between P(z/ a — > i/g) and P{y a — > f a ). The latter can be explicitly 
expressed as follows: 

AP a3 = P( Vp -> i/J - P(i/ a - ^) 

= 16 sin sin F 3 i sin P 32 - sin P 10 sin P 20 sin F 2 i 

- sin Fio sin P 30 sin P 31 ) . (A.34) 

If the hierarchical patterns of neutrino masses and flavor mixing angles are assumed, the 
expression of AP a g may somehow be simplified [143]. Note that only three of the twelve 
nonvanishing asymmetries AP a p are independent, as a consequence of the unitarity of 
V or the correlation of J^. Since only the transition probabilities of active neutrinos 
can be realistically measured, we are more interested in the T-violating asymmetries 
AP e/i , AP^ T and AP re . The overall matter contamination residing in AP a3 is usually 
expected to be insignificant. The reason is simply that the terrestrial matter effects in 
P{y a — > Vq) and P{yg — * v a ), which both depend on the parameters (a, a'), may partly 
(even essentially) cancel each other in the T-violating asymmetry AP^. In contrast, 
P(u a — > u 8 ) and P(v a — > V : /3 ) are associated respectively with (+a, +a') and (—a, —a'), 
thus there should not have large cancellation of matter effects in the corresponding CP- 
violating asymmetries. 



References 

[1] S. Weinberg, Phys. Rev. Lett. 19, 1264 (1967); A. Salam, in Elementary Particle 
Theory, edited by N. Svartholm (Almquist and Wiksells, Stockholm, 1969), p. 367; 
S.L. Glashow, J. Iliopoulos, and L. Maiani, Phys. Rev. D 2, 1285 (1970). 

[2] Y. Fukuda et al, Phys. Lett. B 436, 33 (1998); Phys. Rev. Lett. 81, 1562 (1998). 
For a review, see: C.K. Jung, C. McGrew, T. Kajita, and T. Mann, Ann. Rev. Nucl. 
Part. Sci. 51, 451 (2001). 

[3] SNO Collaboration, Q.R. Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002); Phys. 
Rev. Lett. 89, 011302 (2002). 

[4] KamLAND Collaboration, K. Eguchi et al, Phys. Rev. Lett. 90, 021802 (2003). 

19 Note that the differences of effective neutrino masses A^p (for i = 1,2,3), which must be CP- 
conserving, keep unchanged for the replacements =^> — J^,. Therefore the sign flip of results 

in that of J^g, as indicated by the sum rules given in (A. 32). 



67 



[5] K2K Collaboration, M.H. Aim et al, Phys. Rev. Lett. 90, 041801 (2003). 

[6] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 

[7] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B 420, 397 (1998); Palo 
Verde Collaboration, F. Boehm et al, Phys. Rev. Lett. 84, 3764 (2000). 

[8] LSND Collaboration, C. Athanassopoulos, et al, Phys. Rev. Lett. 81, 1774 (1998); 
Phys. Rev. C 58, 2489 (1998). 

[9] E. Witten, Nucl. Phys. B (Proc. Suppl.) 91, 3 (2001). 

[10] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. 
Theor. Phys. 49, 652 (1973). 

[11] E. Majorana, Nuovo Cim. 14, 171 (1937). 

[12] H. Fritzsch and Z.Z. Xing, Prog. Part. Nucl. Phys. 45, 1 (2000); hep-ph/9912358. 

[13] S. Weinberg, Phys. Rev. D 5, 1962 (1972); A. Zee, Phys, Lett. B 93, 389 (1980). 

[14] H. Fritzsch, Phys. Lett. B 70, 436 (1977); Phys. Lett. B 73, 317 (1978); S. Weinberg, 
in Transactions of the New York Academy of Sciences, 38, 185 (1977); F. Wilczek 
and A. Zee, Phys. Lett. B 70, 418 (1977). 

[15] CD. Froggatt and H.B. Nielsen, Nucl. Phys. B 147, 277 (1979); H. Harari, H. Haut, 
and J. Weyers, Phys. Lett. B 78, 459 (1978). 

[16] T. Yanagida, in Proceedings of the Workshop on Unified Theory and the Baryon 
Number of the Universe, edited by O. Sawada and A. Sugamoto (KEK, 1979); M. 
Gell-Mann, P. Ramond, and R. Slansky, in Super gravity, edited by F. van Nieuwen- 
huizen and D. Freedman (North Holland, Amsterdam, 1979); R.N. Mohapatra and 
G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[17] K.R. Dienes, E. Dudas, and T. Gherghetta, Nucl. Phys. B 557, 25 (1999); N. 
Arkani-Hamed, S. Dimopoulos, G. Dvali, and J. March-Russell, hep-ph/9811448; 
G. Dvali and A.Yu. Smirnov, Nucl. Phys. B 563, 63 (1999). 

[18] G. Altarelli and F. Feruglio, Phys. Rept. 320, 295(1999); hep-ph/0206077; S.M. 
Barr and I. Dorsner, Nucl. Phys. B 585, 79 (2000); S.F. King, hep-ph/0208266. 

[19] R.N. Mohapatra, lectures delivered at the NEUPAST school, Abdus Salam ICTP, 
Trieste, October 2002; hep-ph/0211252; Y. Grossman, lectures given at the TASI 
summer school, University of Colorado, June 2002; hep-ph/0305245. 

[20] A. Lukas, P. Ramond, A. Romanino, and G.G. Ross, JHEP 0104, 010 (2001); and 
references therein. 

[21] H. Murayama, Int. J. Mod. Phys. A 17, 3403 (2002); M.C. Gonzalez-Garcia and Y, 
Nir, Rev. Mod. Phys. 75, 345 (2003); G.G. Ross and L. Velasco-Sevilla, Nucl. Phys. 
B 653, 3 (2003); W. Grimus, hep-ph/0307149. 



68 



[22] H. Fritzsch and P. Minkowski, Ann. Phys. 93, 193 (1975); H. Georgi, in Particles 
and Fields (AIP, New York, 1975). 

[23] J.A. Casas and A. Ibarra, Nucl. Phys. B 618, 171 (2002). 

[24] J. Ellis and M. Raidal, Nucl. Phys. B 643, 229 (2002). 

[25] R.N. Mohapatra and P.B. Pal, Massive Neutrinos in Physics and Astrophysics, 
second edition (World Scientific, 1998). 

[26] PH. Frampton, S.L. Glashow, and T. Yanagida, Phys. Lett. B 548, 119 (2002). 

[27] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

[28] G.L. Fogli, E. Lisi, A. Marrone, and D. Montanino, Phys. Rev. D 67, 093006 (2003). 

[29] Homestake Collaboration, B.T. Cleveland et al, Astrophys. J. 496, 505 (1998). 

[30] GALLEX Collaboration, W. Hampel et al, Phys. Lett. B 447, 127 (1999); GNO 
Collaboration, M. Altmann et al, Phys. Lett. B 490, 16 (2000). 

[31] SAGE Collaboration, J.N. Abdurashitov et al, J. Exp. Theor. Phys. 95, 181 (2002). 

[32] J.N. Bahcall, M.H. Pinsonneault, and S. Basu, Astrophys. J. 555, 990 (2001). 

[33] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); S.P Mikheyev and AYu. Smirnov, 
Sov. J. Nucl. Phys. 42, 913 (1985). For a review, see: T.K. Kuo and J. Pantaleone, 
Rev. Mod. Phys. 61, 937 (1989). 

[34] See, e.g., J.N. Bahcall, P.I. Krastev, and A.Yu. Smirnov, JHEP 05, 015 (2001); and 
references therein. 

[35] V. Barger and D. Marfatia, Phys. Lett. B 555, 144 (2003); G.L. Fogli et al, Phys. 
Rev. D 67, 073002 (2003); M. Maltoni, T. Schwetz, and J.W.F. Valle, Phys. Rev. D 
67, 093003 (2003); J.N. Bahcall, M.C. Gonzalez- Garcia, and C. Pena-Garay, JHEP 
0302, 009 (2003); PC. de Holanda and A.Yu. Smirnov, JCAP 0302, 001 (2003). 

[36] Z.Z. Xing, Phys. Rev. D 65, 077302 (2002). 

[37] M. Aoki et al, Phys. Rev. D 67, 093004 (2003). 

[38] KARMEN Collaboration, K. Eitel et al, Nucl. Phys. B (Proc. Suppl.) 77, 212 
(1999); 87, 281 (2000). 

[39] MiniBooNE Collaboration, R. Tayloe, talk given at Neutrino 2002, Munich, 2002. 

[40] A.D. Dolgov, Phys. Rept. 370, 333 (2002). 

[41] K. Kainulainen and K.A. Olive, hep-ph/0206163; G.G. Raffelt, astro-ph/0207220; 
S.F. King, hep-ph/0210089. 

[42] H. Pas and T.J. Weiler, Phys. Rev. D 63, 113015 (2001); H. Pas, hep-ph/0209018. 



69 



[43] S.M. Bilenky, C. Giunti, J.A. Grifols, and E. Masso, Phys. Rept. 379, 69 (2003). 

[44] C.L. Bennett et al, astro-ph/0302207; D.N. Spergel et al, astro-ph/0302209. 

[45] S.S. Gerstein and Y.B. Zeldovich, JETP Lett. 4, 120 (1966); R. Cowsik and J. 
McClelland, Phys. Rev. Lett. 29, 669 (1972). 

[46] M. Colles et al, MNRAS 328, 1039 (2001). 

[47] S. Hannestad, astro-ph/0303076. 

[48] C. Giunti, Mod. Phys. Lett. A 18, 1179 (2003); A.S. Joshipura and S. Mohanty, 
Phys. Rev. D 67, 091302 (2003); K. Cheung, Phys. Lett. B 562, 97 (2003); G. 
Bhattacharyya, H. Pas, L. Song, and T.J. Weiler, Phys. Lett. B 564, 175 (2003); 
O. Elgaroy and O. Lahav, JCAP 0304, 004 (2003); A. Pierce and H. Murayama, 
hep-ph/0302131. 

[49] Z.Z. Xing, hep-ph/0303178; to appear in HEP&NP. 

[50] KATRIN Collaboration, A. Osipowicz et al, hep-ex/0109033. 

[51] O. Cremonesi, hep-ex/0210007; talk given at the XXth International Conference on 
Neutrino Physics and Astrophysics, May 2002, Munich, Germany. 

[52] H.V. Klapdor-Kleingrothaus, Nucl. Phys. Proc. Suppl. 100, 309 (2001); and refer- 
ences therein. 

[53] C.E. Aalseth et al, Phys. Atom. Nucl. 63, 1268 (2000) [Yad. Fiz. 63, 1341 (2000)]. 

[54] H.V. Klapdor-Kleingrothaus, A. Dietz, H.L. Harney, and I.V. Krivosheina, Mod. 
Phys. Lett. A 16, 2409 (2002). 

[55] F. Feruglio, A. Strumia, and F. Vissani, Nucl. Phys. B 637, 345 (2002); C.E. Aalseth 
et al, Mod. Phys. Lett. A 17, 1475 (2002); Y.G. Zdesenko, F.A. Danevich, and V.I. 
Tretyak, Phys. Lett. B 546, 206 (2002). 

[56] H.V. Klapdor-Kleingrothaus, hep-ph/0302237; hep-ph/0303217; and references 
therein. 

[57] L. Wolfenstein, in Proc. of Neutrino 84, p. 730; F. Vissani, JHEP 9906, 022 (1999); 
S.M. Bilenky, S. Pascoli, and S.T. Petcov, Phys. Rev. D 64, 053010 (2001); W. 
Rodejohann, Nucl. Phys. B 597, 110 (2001); F. Feruglio, A. Strumia, and F. Vissani, 
in Ref. [55]; S. Pascoli, S.T. Petcov, and L. Wolfenstein, Phys. Lett. B 524, 319 
(2002). 

[58] Z.Z. Xing, hep-ph/0305195; to appear in Phys. Rev. D. 

[59] S. Pascoli and S.T. Petcov, Phys. Lett. B 544, 239 (2002); S. Pascoli, S.T. Petcov, 
and W. Rodejohann, Phys. Lett. B 549, 177 (2002); F.R. Joaquim, hep-ph/0304276. 

[60] PH. Frampton, S.L. Glashow, and D. Marfatia, Phys. Lett. B 536, 79 (2002). 



70 



[61] H. Fritzsch and Z.Z. Xing, Phys. Lett. B 372, 265 (1996); Phys. Lett. B 440, 313 

(1998) ; Phys. Rev. D 61, 073016 (2000). 

[62] We thank W.L. Guo for a careful analysis of all twenty three-zero textures of the 
neutrino mass matrix. A similar conclusion has also been drawn in Ref. [60]. 

[63] Z.Z. Xing, Phys. Lett. B 530, 159 (2002); Phys. Lett. B 539, 85 (2002); W.L. Guo 
and Z.Z. Xing, Phys. Rev. D 67, 053002 (2003). 

[64] Z.Z. Xing, hep-ph/0307007. 

[65] H. Fritzsch, Nucl. Phys. B 155, 189 (1979). 

[66] Z.Z. Xing, Phys. Lett. B 550, 178 (2002). 

[67] Z.Z. Xing and H. Zhang, hep-ph/0304234; to appear in Phys. Lett. B. 

[68] See, e.g., D. Du and Z.Z. Xing, Phys. Rev. D 48, 2349 (1993); H. Fritzsch and 
Z.Z. Xing, Phys. Lett. B 353, 114 (1995); Nucl. Phys. B 556, 49 (1999); Phys. 
Lett. B 506, 109 (2001); Phys. Lett. B 555, 63 (2003); PS. Gill and M. Gupta, 
Phys. Rev. D 56, 3143 (1997); Phys. Rev. D 57, 3971 (1998); G.C. Branco, D. 
Emmanuel-Costa, and R.G. Felipe, Phys. Lett. B 477, 147 (2000); R. Rosenfeld 
and J.L. Rosner, Phys. Lett. B 516, 408 (2001); J.L. Chkareuli and CD. Froggatt, 
Nucl. Phys. B 626, 307 (2002); JW. Mei and Z.Z. Xing, Phys. Rev. D 67, 077301 
(2003); and references therein. 

[69] For recent reviews with extensive references, see: S. Raby, hep-ph/9501349 (unpub- 
lished); M.C. Chen and K.T. Mahanthappa, hep-ph/0305088; G. Altarelli and F. 
Feruglio, hep-ph/0306265. 

[70] See, e.g., H. Nishiura, K. Matsuda, and T. Fukuyama, Phys. Rev. D 60, 013006 

(1999) ; J.L. Chkareuli and CD. Froggatt, Phys. Lett. B 450, 158 (1999); H. Fritzsch 
and Z.Z. Xing, in Ref. [12]; K. Matsuda, T. Fukuyama, and H. Nishiura, Phys. Rev. 
D 61, 053001 (2000); W. Buchmuller and D. Wyler, Phys. Lett. B 521, 291 (2001); 
Z.Z. Xing, Phys. Lett. B 545, 352 (2002); hep-ph/0209066; M. Bando and M. Obara, 
hep-ph/0212242; Prog. Theor. Phys. 109, 995 (2003); D. Falcone, hep-ph/0303074. 

[71] A. Kageyama, S. Kaneko, N. Shimoyama, and M. Tanimoto, Phys. Lett. B 538, 
96 (2002); R. Barbieri, T. Hambye, and A. Romanino, JHEP 0303, 017 (2003); M. 
Honda, S. Kaneko, and M. Tanimoto, hep-ph/0303227. 

[72] S.M. Bilenky, J. Hosek, and S.T. Petcov, Phys. Lett. B 94, 495 (1980); J. Schechter 
and J.W.F. Valle, Phys. Rev. D 22, 2227 (1980). 

[73] Z.Z. Xing, J. Phys. G 28, B7 (2002); hep-ph/0112120. 

[74] M. Fukugita and M. Tanimoto, Phys. Lett. B 515, 30 (2001). 

[75] L. Maiani, in Proceedings of the International Symposium on lepton and Photon 
Interactions at High Energies (DESY, Hamburg, 1977), p. 867; L.L. Chau and 
W.Y. Keung, Phys. Rev. Lett. 53, 1802 (1984); H. Fritzsch, Phys. Rev. D 32, 3058 
(1985); H. Harari and M. Leurer, Phys. Lett. B 181, 123 (1986). 



71 



[76] W.L. Guo and Z.Z. Xing, Phys. Rev. D 67, 053002 (2003). 
[77] F.J. Botella and L.L. Chau, Phys. Lett. B 168, 97 (1986). 

[78] H. Fritzsch and Z.Z. Xing, Phys. Lett. B 413, 396 (1997); Phys. Rev. D 57, 594 

(1998) . 

[79] H. Fritzsch and Z.Z. Xing, Phys. Lett. B 517, 363 (2001). 

[80] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985); D.D. Wu, Phys. Rev. D 33, 860 
(1986). 

[81] A. de Gouvea, B. Kayser, and R.N. Mohapatra, Phys. Rev. D 67, 053004 (2003). 
[82] N. Cabibbo, Phys. Lett. B 72, 333 (1978). 

[83] See, e.g., M. Tanimoto, T. Watari, and T. Yanagida, Phys. Lett. B 461, 345 (1999); 
M. Tanimoto, Phys. Lett. B 483, 417 (2000); Y. Koide and A. Ghosal, Phys. Lett. B 
488, 344 (2000); T. Miura, T. Shindou, E. Takasugi, and M. Yoshimura, Phys. Rev. 
D 63, 053006 (2001); Z.Z. Xing, Phys. Rev. D 63, 057301 (2001); E.K. Akhmedov, 
G.C. Branco, F.R. Joaquim, and J.I. Silva-Marcos, Phys. Lett. B 498, 237 (2001); 
G.C. Branco and J.I. Silva-Marcos, Phys. Lett. B 526, 104 (2002); J. Kubo, A. 
Mondragon, M. Mondragon, and E. Rodriguez- Jauregui, Prog. Theor. Phys. 109, 
795 (2003). 

[84] V. Barger, S. Pakvasa, T.J. Weiler, and K. Whisnant, Phys. Lett. B 437, 107 (1998). 
[85] F. Vissani, hep-ph/9708483 (unpublished). 

[86] See, e.g., M. Jezabek and Y. Sumino, Phys. Lett. B 440, 327 (1998); E. Ma, D.P 
Roy, and U. Sarkar, Phys. Lett. B 444, 391 (1998); S. Davidson and S.F. King, Phys. 
Lett. B 445, 191 (1998); Y. Nomura and T. Yanagida, Phys. Rev. D 59, 017303 

(1999) ; C. Jarlskog, M. Matsuda, S. Skadhauge, and M. Tanimoto, Phys. Lett. B 
449, 240 (1999); C.H. Albright and S.M. Barr, Phys. Lett. B 461, 218 (1999); S.K. 
Kang and C.S. Kim, Phys. Rev. D 59, 091302 (1999); Y.L. Wu, Eur. Phys. J. C 10, 
491 (1999); Sci. China A 43, 988 (2000); Q. Shafi and Z. Tavartkiladze, Phys. Lett. 
B 482, 145 (2000); Phys. Lett. B 487, 145 (2000); C.S. Kim and J.D. Kim, Phys. 
Rev. D 61, 057302 (2000); Z.Z. Xing, Phys. Rev. D 61, 057301 (2000); Phys. Rev. 
D 64, 017304 (2001); Phys. Rev. D 64, 093013 (2001); D.V. Ahluwalia, Mod. Phys. 
Lett. A 16, 917 (2001); Y. Koide and A. Ghosal, Phys. Rev. D 63, 037301 (2001); 
W. Grimus and L. Lavoura, JHEP 0107, 045 (2001); hep-ph/0305046; C.S. Lam, 
Phys. Lett. B 507, 214 (2002); T. Kitabayashi and M. Yasue, Phys. Lett. B 524, 
308 (2002); Int. J. Mod. Phys. A 17, 2519 (2002); K.S. Babu and S.M. Barr, Phys. 
Lett. B 525, 289 (2002); S. Antusch, J. Kersten, M. Lindner, and M. Ratz, Phys. 
Lett. B 544, 1 (2002); C. Giunti and M. Tanimoto, Phys. Rev. D 66, 113006 (2002); 
T. Ohlsson and G. Seidl, Nucl. Phys. B 643, 247 (2002); M.C. Gonzalez- Garcia, 
hep-ph/0210359. 

[87] L. Wolfenstein, Phys. Rev. D 18, 958 (1978). 



72 



[88] P.F. Harrison, D.H. Perkins, and W.G. Scott, Phys. Lett. B 530, 167 (2002); Z.Z. 
Xing, Phys. Lett. B 533, 85 (2002). X.G. He and A. Zee, Phys. Lett. B 560, 87 
(2003); C.I. Low and R.R. Volkas, hep-ph/0305243; A. Zee, hep-ph/0307323. 

[89] C. Giunti, hep-ph/0209103. 

[90] R. Peccei, lectures given at Topical Seminars on Frontier of Particle Physics: Neu- 
trinos and Cosmology, August 2002, Beijing. 

[91] Z.Z. Xing, hep-ph/0211465; to appear in J. Phys. G. 

[92] For recent reviews with extensive references, see: J. A. Casas, J.R. Espinosa, A. 
Ibarra, and I. Navarro, Nucl. Phys. B 573, 652 (2000); PH. Chankowski and S. 
Pokorski, Int. J. Mod. Phys. A 17, 575 (2002); S. Antusch, J. Kersten, M. Lindner, 
and M. Ratz, hep-ph/0305273. 

[93] Z.Z. Xing, Phys. Rev. D 63, 073012 (2001). 

[94] I. Mocioiu and R. Shrock, Phys. Rev. D 62, 053017 (2000); and references therein. 

[95] V. Barger, K. Whisnant, S. Pakvasa, and R.J.N. Phillips, Phys. Rev. D 22, 2718 
(1980); H.W. Zaglauer and K.H. Schwarzer, Z. Phys. C 40, 273 (1988). 

[96] J. Pantaleone, Phys. Rev. D 49, 2152 (1994); V. Barger, S. Geer, R. Raja, and K. 
Whisnant, Phys. Rev. D 62, 013004 (2000). 

[97] Z.Z. Xing, Phys. Lett. B 487, 327 (2000). 

[98] Z.Z. Xing, Phys. Rev. D 64, 033005 (2001); Fortsch. Phys. 50, 569 (2002). 

[99] K. Kimura, A. Takamura, and H. Yokomakura, Phys. Lett. B 537, 86 (2002); Phys. 
Rev. D 66, 073005 (2002). 

[100] T. Ohlsson and H. Snellman, J. Math. Phys. 41, 2768 (2000); Phys. Lett. B 474, 
153 (2000); P.F. Harrison and W.G. Scott, Phys. Lett. B 535, 229 (2002); P.F. 
Harrison, W.G. Scott, and T.J. Weiler, Phys. Lett. B 565, 159 (2003). 

[101] Z.Z. Xing, Phys. Rev. D 64, 073014 (2001). 

[102] S. Toshev, Mod. Phys. Lett. A 6, 455 (1991). 

[103] M. Freund, Phys. Rev. D 64, 053003 (2001). 

[104] A. Dziewonski, Earth Structure, in "The Encyclopedia of Solid Earch Geophysics", 
edited by D.E. James (Van Nostrand Reinhold, New York, 1989), p. 331; I. Mocioiu 
and R. Shrock, in Ref. [94]. 

[105] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

[106] Modified versions of the Wolfenstein parametrization with higher accuracy can be 
found, e.g., in: A.J. Buras, M.E. Lautenbacher, and G. Ostermaier, Phys. Rev. D 
50, 3433 (1994); Z.Z. Xing, Phys. Rev. D 51, 3958 (1995); P. Paganini, F. Parodi, 
P. Roudeau, and A. Stocchi, Physica Scripta 58, 556 (1998). 



73 



Z.Z. Xing, Phys. Rev. D 65, 113010 (2002). 

Z.Z. Xing, Nucl. Phys. B (Proc. Suppl.) 50, 24 (1996); Nuovo Cim. A 109, 115 
(1996); J. Phys. G 23, 717 (1997). 

P. Kaus and S. Meshkov, hep-ph/0211338. 

K. Sasaki, Z. Phys. C 32, 149 (1986). 

V.A. Naumov, Int. J. Mod. Phys. D 1, 379 (1992). 

PF. Harrison and W.G. Scott, Phys. Lett. B 476, 349 (2000). 

H. Fritzsch and Z.Z. Xing, Nucl. Phys. B 556, 49 (1999). 

H.S. Chen et al, hep-ph/0 104266; and references therein. 

See, e.g., J. Sato, Nucl. Instrum. Meth. A 472, 434 (2000). 

Y. Farzan and A.Yu. Smirnov, Phys. Rev. D 65, 113001 (2002). 

A.D. Sakharov, JETP Lett. 5, 24 (1967). 

For a recent review with extensive references, see: M. Dine and A. Kusenko, hep- 
ph/0303065 (to appear in Rev. Mod. Phys.). 

M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 

See, e.g., S. Davidson and A. Ibarra, Phys. Lett. B 535, 25 (2002). 

M.A. Luty, Phys. Rev. D 45, 455 (1992); M. Flanz, E.A. Paschos, and U. Sarkar, 
Phys. Lett. B 345, 248 (1995); L. Covi, E. Roulet, and F. Vissani, Phys. Lett. B 
384, 169 (1996); M. Flanz, E.A. Paschos, U. Sarkar, and J. Wess, Phys. Lett. B 
389, 693 (1996); M. Pliimacher, Z. Phys. C 74, 549 (1997) A. Pilaftsis, Phys. Rev. 
D 56, 5431 (1997); Int. J. Mod. Phys. A 14, 1811 (1999); W. Buchmuller and M. 
Pliimacher, Phys. Lett. B 431, 354 (1998); M. Pliimacher, Nucl. Phys. B 530, 207 
(1998); R. Barbieri, P. Creminelli, A. Strumia, and N. Tetradis, Nucl. Phys. B 575, 
61 (2000); T. Hambye, Nucl. Phys. B 633, 171 (2002); W. Buchmuller, P. Di Bari, 
and M. Pliimacher, Nucl. Phys. B 643, 367 (2002). 

[122] G. 't Hooft, Phys. Rev. Lett. 37, 8 (1976); V.A. Kuzmin, V.A. Rubakov, and M.E. 
Shaposhnikov, Phys. Lett. B 155, 36 (1985). 

[123] J.A. Harvey and M.S. Turner, Phys. Rev. D 42, 3344 (1990); SYu. Khlebnikov 
and M.E. Shaposhnikov, Nucl. Phys. B 308, 169 (1988). See also, K. Hamaguchi, 
PhD thesis, hep-ph/0212305. 

[124] E.W. Kolb and M.S. Turner, The Early Universe, Addison- Wesley (1990); M. 
Pliimacher, Ph.D. Thesis, hep-ph/9807557. 

[125] E.W. Kolb and M.S. Turner, in Ref. [124]; H.B. Nielsen and Y. Takanishi, Phys. 
Lett. B 507, 241 (2001); M. Hirsch and S.F. King, Phys. Rev. D 64, 113005 (2001). 

74 



[107 
[108 

[109 
[110 
[111 
[112; 
[113 
[114; 
[115 
[116 
[117; 
[118 

[119 
[120 
[12k 



[126] J. Ellis, M. Raidal, and T. Yanagida, Phys. Lett. B 546, 228 (2002); and references 
therein. 

[127] See, e.g., G.F. Giudice, M. Peloso, A. Riotto, and I. Tkachev, JHEP 9908, 014 
(1999); T. Asaka, K. Hamaguchi, M. Kawasaki, and T. Yanagida, Phys. Lett. B 
464, 12 (1999); Phys. Rev. D 61, 083512 (2000); J. Garcia-Bellido and E. Ruiz 
Morales, Phys. Lett. B 536, 193 (2002); M. Fujii, K. Hamaguchi, and T. Yanagida, 
Phys. Rev. D 65, 115012 (2002). 

[128] Z.Z. Xing, Phys. Lett. B 545, 352 (2002). 

[129] W. Buchmiiller and M. Pliimacher, Phys. Lett. B 389, 73 (1996). 

[130] Z.Z. Xing, in Proceedings of ICHEP 2002, edited by S. Bentvelsen, P. de Jong, J. 
Koch, and E. Laenen (Elsevier Science, 2003), p. 43; hep-ph/0209066. 

[131] See, e.g., G.C. Branco, T. Morozumi, B.M. Nobre, and M.N. Rebelo, Nucl. Phys. 
B 617, 475 (2001); D. Falcone and F. Tramontano, Phys. Rev. D 63, 073007 (2001); 
G.C. Branco, R. Gonzalez Felipe, F.R. Joaquim, and M.N. Rebelo, Nucl. Phys. B 
640, 202 (2002); M.N. Rebelo, Phys. Rev. D 67, 013008 (2003); S. Pascoli, S.T. 
Petcov, and W. Rodejohann, hep-ph/0302054; W. Buchmiiller, P. Di Bari, and M. 
Pliimacher, hep-ph/0302092; D. Falcone, hep-ph/0305229. 

[132] W. Buchmiiller and D. Wyler, Phys. Lett. B 521, 291 (2001); PH. Frampton, 
S.L. Glashow, and T. Yanagida, in Ref. [26]; T. Endoh, S. Kaneko, S.K. Kang, T. 
Morozumi, and M. Tanimoto, Phys. Rev. Lett. 89, 231601 (2002); S.F. King, Phys. 
Rev. D 67, 113010 (2003). 

[133] W. Buchmiiller and M. Pliimacher, Int. J. Mod. Phys. A 15, 5047 (2000); W. 
Buchmiiller, hep-ph/0204288; J.C. Pati, hep-ph/0209160. 

[134] H. Murayama, hep-ph/0208005. 

[135] H. Murayama and T. Yanagida, Phys. Lett. B 520, 263 (2001); G. Barenboim, L. 
Borissov, and J. Lykken, hep-ph/0212116. 

[136] See, e.g., K.S. Babu and S. Pakavasa, hep-ph/0204236. 

[137] A. Strumia, Phys. Lett. B 539, 91 (2002); M. Maltoni, T. Schwetz, M.A. Tortola, 
and J.W.F. Valle, Nucl. Phys. B 643, 321 (2002); Nucl. Phys. Proc. Suppl. 114, 
203 (2003); PC. de Holanda and A.Yu. Smirnov, hep-ph/0307266. 

[138] V. Barger, Y.B. Dai, K. Whisnant, and B.L. Young, Phys. Rev. D 59, 113010 
(1999). 

[139] H. Fritzsch and J. Plankl, Phys. Rev. D 35, 1732 (1987); Z.Z. Xing and D.Du, 
Chin. Phys. Lett. 9, 9 (1992). 

[140] W.L. Guo and Z.Z. Xing, Phys. Rev. D 65, 073020 (2002). 

[141] W.L. Guo and Z.Z. Xing, Phys. Rev. D 66, 097302 (2002). 



75 



[142] S.M. Bilenky, C. Giunti, and W. Grimus, Prog. Part. Nucl. Phys. 43, 1 (1999). 
[143] A. Kalliomaki, J. Maalampi, and M. Tanimoto, Phys. Lett. B 469, 179 (1999). 



76 



